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I.  FUNDAMENTAL  AND  DERIVED  UNITS 


MctHe 

'  ' 

. . . . — . . . 

English 

Symbol 

Unit 

Abhrevla- 

tlon 

Unit 

Ahtirevla- 

tion 

length »  - _ 

r 

meter . . . 

ft,  (or  ml.) 
sec.  (or  hr.) 

Time.. _ 

t 

second.  . . .  . . 

A 

second  (or  hour) . 

Korea.. . 

F 

weight  of  1  kilogram..... 

kg 

weight  of  1  (Keind .... 

lb. 

I'ower.. _ _ _ 

P 

horaopower  (metrie) ..... 

horsepower,...  . . 

hji. 

, .... 

v 

/kilometers  per  hour...... 

k.pb. 

miles  per  hour. ... 

m.p.h. 

(meters  per  aeoond . 

mp.a, 

fret  |»er  second...  . 

f.p.s. 

2.  GENERAL  SYMBOLS 


W,  Weight---  mfl 

q ,  Standard  acceleration  of  gravity  «>0.80(I65 

in/*1  or  32.1740  ft./Hoc.* 
u' 


/,  Moment  of  inorliarrmi'.  (Indicate  axis  of 
mdiirn  of  (gyration  k  by  proper  Bulwcript.) 

( ’oolficiont  of  viacoaily 


v,  Kinematic  viscosity 

Density  (mtM  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m ‘-s*  ut 
15°  C.  and  700  urn;;  or  0.002378  Ib.-ft  ace.* 
Specific  weight  of  “standard”  air,  1.2255  kg/tn1  or 
0.07651  ll»./cu.  ft. 


3.  AERODYNAMIC  SYMBOLS 


", 

■S’., 

O, 

f>, 

f, 

lj‘ 

s’ 

V, 

V , 


is, 

is, 

V, 


n, 


Area 

Area  of  wing 
( inp 
Spun 
(  tionl 


/><  i»  ‘  I  ratio 
True  idr  speed 
I  lynumic  prcaauro^’^pT7* 

I, (ft,  absolute  coefficient 

Drag,  absolute  coeffiejent 

Profile  drag,  almoluto  coefficient  Cn,** 

Induced  drag,  absolute  coefficient  Cd,1 

far, mi  to  drag,  abfloluto  coefficient 

’  r/o 

I  ro^  wind  force,  absolute  coefficient  Cn^^L 


3 


Q‘ 


T) 


3 


t»,  Anglo  of  sotting  of  wings  (relative  to  thrust 
line) 

i„  Anglo  of  stabilizer  setting  (relative  to  thrust 

line) 

Q,  Rcsuhant  moment 

O,  Resultant  angular  velocity 

Vl 

p  — ,  Reynold#  Number,  whore  fi#  a  linear  dimension 
M  (c.g..  for  a  model  airfoil  3  In.  chord,  100 
m.p.n.  normal  pramuro  at  15*  C.,  tho  eor- 
r  or  ponding  numlmr  is  234,000;  or  for  a  model 
of  10  cm  chord,  40  m.p.a.,  tho  corresponding 
number  la  274,000) 

Cf,  Ccnfer-of-presaure  coefficient  (ratio  of  distance 

of  c.p.  from  leading  edge  to  chord  length) 
a,  Angie  of  attack 

«,  Angle  of  dowmraalt 

A#,  Angle  of  attack,  infinite  aspect  ratio 

n„  Angle  of  attack,  induced 

«„  Angle  of  attack,  absolute  tmmsumd  from  mo- 
lift  position) 
y,  Flight-path  angle 
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MECHANISM  OF  FLUTTER 

A  THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION  OF  THE 

FLUTTER  PROBLEM 

By  Thkodokk  Tiikiidoksen  and  1.  K.  Cakhii  k 


SUMMARY 

The  results  of  the  basic  flutter  theory  originally  devised 
in  ]!)34  anil  published  ax  X.  .1.  t *1.  Technical  Report 
Xo.  .{.W  are  presented  in  a  simpler  and  more  complete 
form  con  celt  lent  for  further  studies.  The  paper  attempts 
to  facilitate  the  judgment  of  flutter  problems  by  a  syste¬ 
matic  xiii rey  of  the  theoretical  effects  oj  the  carious  param¬ 
eters.  A  large  number  of  experiments  were  conducted  on 
ca  utile  re  r  wings,  with  anil  without  ailerons,  in  the 
A.  .1.  ('.  .1.  high-speed  wind  tunnel  for  the  purpose  of 
ceri.fying  the  theory  and  to  study  its  adaptability  to  three- 
dimensional  problems.  The  experiments  included  studies 
on  wing  taper  ratios,  nacelles,  attached  floats,  and  external 
bracings.  The  essential  effects  in  the  transition  to  the 
three-dimensional  problem  hare  been  established.  Of 
particular  interest  is  the  existence  of  specific  flutter  modes 
as  distinguished  from  ordinary  vibration  miules.  On  the 
basis  of  the  concepts  introeluced,  results  that  are  apjmrently 
paradoxical  coulef  logically  be  brought  into  conformity 
with  the.  theory.  In  fact,  it  is  shown  that  there  exists  a 
r  ithrr  remarkable,  agreement  between  tlieoietical  anil 
experimental  results.  simple  method  is  presented  for 
numerical  calculations  of  the  flutter  speed  by  routine 
operations,  rei/uiring  no  reference  to  the  theory.  Appli¬ 
cation  is  made  to  a  complete  numerical  example.  The 
matter  of  identifying  possible  types  of flutter  in  an  airplane 
ami  of  determining  the  parameters  is  briefly  iliscusseil. 
A  section  treating  the  subject  of  forced  vibrations  of  a 
wing  in  an  ihr  stream  and  the  question  of  air  damping 
in  its  relation  to  flutter  is  included. 

INTRODUCTION 

The  theory  of  flutter.-  The  problem  of  flutter  is 
passing  through  u  period  of  rapid  development.  Full 
cognizance  is  taken  of  the  value  of  the  theory;  a  simple 
or  an  empirical  understanding  of  this  problem  is  not 
available  ami  could,  at  best,  be  of  value  only  to  the 
investigator.  An  exact  treatment  of  the  basic  flutter 
problem  in  two-dimensional  flow,  involving  the  im¬ 
portant  functions  F  and  G  relating  to  the  air  forces, 
was  given  bv  Theodorscn  in  UK14.  (,See  reference  1.) 


These  functions  are  simple  combinations  of  Bessel 
functions;  they  have  been  rederived  in  related  form 
by  Cicala  (reference  2)  in  1  by  Kassner  and 
Ki’.gndo  (reference  :t)  in  ItIUti,  and  also  by  Iviissner 
(reference  4)  in  15M0,  who  pointed  out  the  identity  of 
the  functions.  At  about  this  time,  Garrick  (reference  •'>) 
also  established  a  check  on  the  general  functions  /’and 
G  by  comparing  them  with  expressions  by  Wagner, 
Glauert,  and  von  Kdrmdn  and  Burgers  for  special  cases. 

The  system  of  equations  as  given  in  the  original 
paper  is 

(.1)  o.l  exfJi a(A  Ol )  .  OpF  .  (f^F-  hU 

-2("'  ° 

(«)  all  :  a [l  •  fid  ■  (Sj  K  ■  01.  :  h  .1/+  ^4  J-(U)/'  (I 


<o  o.v  •  «;v> .  or  .  o'aj  ■  hii  ■  hs  lv'Vui/  ii 

o  b  ti 

where  A,  li,  (',  etc.,  are  given  on  page  111  of  reference  1. 


ami  (table  I) 


Butting 


rU-i  t\k- 1  •  it  I  (hi 


u„  < 


when* 


0 

h  /m-'-'v 


the  determinant  of  the  coefficients  of  o0,  (J,„  and  lt„ 
becomes 

\Rna  I  I  Iia  R.,i  ^  t  I  nS  It  ah  llnh 

Ki*}J  lif,a  ♦  1 1  f,a  It  hi  *  tl hfl  Ithh'\  1 1  bh 

It  (-a  !  1 1  ca  It  c, 1 ' 1  1 1  ch~\  1  /  ch 

when*  the  IV s  and  /’s  are  listed  in  tlie  appendix. 
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Tliis  determinant  put  equal  to  zero  contains  two  simul¬ 
taneous  equations,  the  solution  of  which  determines  A' 
and  1  fk  from  which  the  (unknowns)  llutter  velocity  and 
flutter  frequency  are  obtained.  Only  the  diagonal 
terms  with  liars  contain  the  quantity  A'.  All  terms 
contain  or  are  functions  of  1  Ik. 

The  work  of  numerous  investigators,  notably  Becker 
and  Koppl,  has  shown  that  the  structural  friction  is 
mainly  a  function  of  amplitude,  not  of  frequency. 
The  structural  friction  can  he  described  by  a  force  in 
phase  with  the  velocity  but  of  a  magnitude  proportional 
to  the  restoring  force.  With  each  restoring-force 
term,  say  a( there  will  be  a  friction  term  ia  t/a  in 
which  tj„  is  the  damping  coellicicnt.  The  net  result  is 
very  simply  that  the  restoring-force  terms  a( /SCj, 
ami  hi\  have  been  replaced  by  terms  of  the  form  a('a 
0  -f  i<Ja),  fiC’A  1  i  itjs),  A  C»  ( 1  r  <!U)-  These  friction  coef¬ 
ficients  occur  only  in  the  diagonal  I  terms  of  the  deter¬ 
minant. 

Technical  flutter  problems  and  the  flutter  param¬ 
eters. — -Experimental  evidence,  some  of  which  is  pre¬ 
sented  later,  has  been  accumulated  which  indicates 
that,  in  the  two-dimensional  problem,  the  flutter  speed 
can  lie  closely  predicted  from  the  theory  if  (lie  para¬ 
meters  are  given  with  accuracy.  (In  fact,  it  seems  that 
in  some  cases  the  flutter  speed  can  be  used  to  determine 
some  parameters  more  precisely  than  by  a  direct 
method.)  In  the  two-dimensional  problem  of  flexure- 
torsion-aileron  flutter,  about  a  dozen  different  quantities 
are  required  to  calculate  the  flutter  speed.  The 
determination  of  these  parameters  requires  technical 
skill  and  experience  and  is  perhaps  the  most  difficult 
stej)  in  the  solution  of  the  flutter  problem.  A  know  ledge 
of  the  functional  dependency  of  the  flutter  speed  on 
each  parameter  is  essential  in  order  to  obtain  sufficient 
accuracy  in  the  determination  of  the  important  ones 
and  to  prevent  waste  of  time  on  those  of  less  influence. 
This  need  is  partly  the  purpose  of  the  material  given  in 
this  paper. 

One  of  the  problems  in  connection  with  an  actual 
airplane  is  the  identification  of  the  combination  of 
vibration  modes  that  may  cause  flutter.  In  regard 
to  wing  flutter,  in  the  case  of  flexure-torsion,  the  situa¬ 
tion  is  fairly  clear.  It  will  be  shown  that  normally  the 
most  important  paiameter  is  the  center-of-gravity 
location.  This  constant  can  he  obtained  w  ith  consider¬ 
able  accuracy  in  the  design  stage.  An  accurate  value 
of  this  parameter  can  also  be  experimentally  obtained 
as  the  “dynamic’’  torsion  axis,  that  is,  the  axis  around 
which  the  wing,  owing  to  the  low  bending  frequency, 
oscillates  when  put  into  torsional  resonance.  The 
location  of  the  (static)  torsional  stiffness  axis  is  much 
more  difficult  to  calculate  or  to  determine  experimentally 
but  fortunately,  as  will  be  observed,  its  effect  on  the 
flutter  velocity  is  small. 

The  internal  damping  coefficients  are,  moreover,  of 
fairly  small  influence  in  flexure-torsion  flutter;  these 


parameters  are  also  fairly  difficult  to  obtain.  On  the 
whole,  however,  it  may  he  said  that  this  ease  of 
flutter  can  be  fairly  well  handled. 

Another  important  case,  that  of  the  combination 
flexure-aileron,  was  shown  by  the  original  study 
(reference  1)  to  be  an  essentially  different  type  of 
flutter  from  flexure-torsion.  Its  primary  characteristic 
is  that,  normally,  the  flutter  is  limited  to  a  r<nnje  of 
speed.s.  Below  and  above  the  extremes  of  this  speed 
range  there  is  aerodynamic  stability.  A  reduction  of 
the  static  moment  of  the  aileron  with  respect  to  the 
hinge  (balancing)  reduces  the  range;  that  is,  the  lower 
limit  is  raised  and  the  upper  limit  is  lowered.  Damping 
in  the  structure  is  found  to  have  the  same  general  effect. 
Sufficient  internal  friction  will,  in  fact,  completely 
eliminate  the  danger  of  flutter,  as  will  also  complete 
mass  balance  or  the  proper  combination  of  both. 
The  structural  friction  of  a  wing  system,  although  not 
readily  predictable,  can  he  obtained  hv  a  ground  test. 

In  regard  to  the  tail  assembly,  the  difficulty  is  some¬ 
what  greater  since  it  may  not  be  easy  to  identify  the 
most  dangerous  combination  or  to  predict  or  even  to 
measure  the  necessary  parameters,  including  the 
structural  damping.  It  therefore  seems  that  cotaii 
empirical  or  scmicmpirical  aids  will  be  required  ami 
that  it  will  be  necessary,  for  a  time  at  least,  to  resort 
to  flight-test  methods  as  a  final  assurance  against  tail 
flutter. 

The  transition  to  the  three-dimensional  case  of  actual 
flutter  is  quite  complex.  It  is  necessary  to  consider  an 
aerodynamic  span  effect  (which  fortunately  is  very 
small,  see  reference  til:  the  variation  of  the  parameters 
along  the  span;  the  possibility  of  higher-order  deflec¬ 
tion  inodes;  and,  in  certain  cases,  fractional  span  effects, 
as  for  partial  ailerons.  The  most  promising  manner 
of  attack  on  such  problems  is  by  means  of  the  two- 
dimensional  treatment  with  the  introduction  of  certain 
weight  functions  and  average  parameters  in  conjunc¬ 
tion  with  a  study  of  representative  models  of  reason¬ 
able  simplicity,  followed  by  a  crystallization  of  the 
collected  experience  into  generally  applicable  semi- 
empirical  correction  factors.  The  present  paper  makes 
initial  studies  with  this  purpose  in  mind. 

It  is  realized  that,  for  high  values  of  the  flutter  speed, 
a  correction  must  be  made  for  the  effect  of  compressi¬ 
bility.  In  the  first  order,  this  effect  is  due  to  a  change 
in  the  slope  of  the  lift  curve.  The  air  forces  in  the 
steady  case  are  known  to  he  increased  approximately 
in  the  ratio  1  \  1  .!/•',  where  M  is  the  Mach  number. 

Consequently,  a  decrease  in  the  flutter  speed,  roughly 
as  (I  —  .\f*)Ui,  is  expected.  This  correction,  although 
small  through  the  usual  flight  range,  becomes  appreci¬ 
able  for  speeds  near  sound  speed.  Until  experimental 
verification  is  available,  such  correction  is  preferable 
to  none  and  should  bo  applied  for  high-speed  airplanes. 
(See  footnote  2,  p.  9,  for  details.) 
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MECHANISM  OK  FLUTTER 


Content  of  paper.  A  straightforward  scheme  is  pre¬ 
sented  in  the  first  section  for  ii>iitiiie  calculation  of  the 
flutter  speeds  in  the  two-dimensional  types;  ease  1, 
flexure-torsion;  ease  2,  flexure-aileron;  and  ease*  3, 
torsion-aileron;  and,  in  the  three-degrccs-of-freedom 
type,  flexure-torsion-aileron.  A  numerical  example, 
referring  to  a  modem  large  airplane,  is  included. 

The  second  section  deals  with  a  survey  of  the  effect 
of  the  flutter  parameters  on  the  critical  velocity. 
The  effect  of  changing  the  parameters  within  certain 
practical  limits  in  eases  1, 2,  and  3  is  shown  by  a  num¬ 
ber  of  charts. 

The  discussion  in  the  third  section  deals  with  the 
transition  to  a  three-dimensional  case,  showing  how  a 
"representative”  two-dimensional  wing  may  he  used 
to  (live  the  essential  results.  Both  uniform  and  t altered 
cantilever  wings  are  included.  The  question  of  the 
probable  occurrence  of  higher-order  bonding  modes  in 
flutter  is  also  discussed.  The  effects  of  "friction"  and 
“coupling"  are  especially  pronounced  in  higher-order 
flutter. 

It  is  [minted  out  that  the  deflection  mode  occurring 
in  flutter  is  quite  different  from  that  of  the  static 
condition  and  that  the  lowest  bending  frequency 
involved  in  flutter  is  greater  than  that  of  the  lowest 
ordinary  vibration  mode.  A  new  concept  of  flutter, 
that  the  mode  arising  in  flutter  is  such  that  the  flutter 
speed  is  a  minimum,  is  then  introduced.  In  other 
words,  if  nl!  primary  variables  including  friction  could 
be  included  in  the  analysis,  the  actual  inode  would  be 
determined  from  all  possible  inodes  as  the  one  giving 
the  minimum  critical  speed.  This  concept  is  useful  in 
explaining  certain  otherwise  paradoxical,  experimental 
facts.  The  extreme  difficulty  of  a  direct  analytic 
attack  on  the  general  case,  even  if  all  the  physical 
parameters  were  specified,  justifies  the  adaptin':  of  the 
two-dimensional  treatment  supplemented  by  empirical 
information  obtained  on  actual  wings.  In  fact,  as  will 
later  be  shown,  the  corrections  are  small. 

Almost  100  separate  experiments  were  conducted 
in  the  8-foot  high-speed  tunnel.  The  fourth  section 
deals  with  the  experimental  tests  and  results.  About 
one-half  of  these  tests  pertain  to  flutter  of  wings  in 
flexure-torsion;  the  rest  pertain  to  aileron  flutter. 
Cantilever  wings  of  aluminum  and  of  built-up  wood 
construction  were  used.  The  tests  were  performed  on 
a  conveniently  large  scale,  most  of  the  wings  having  a 
chord  of  I  font  and  a  span  of  a Immi t  7  feet.  The  air 
speeds  ranged  from  .*>()  to  iilauit  300  miles  jht  hour- 
A  number  of  safety  devices  had  to  be  employed  to 
prevent  the  ruin  of  the  tunnel  equipment. 

A  section  is  included  showing  the  theoretical  effects 
of  the  air  damping  on  the  forced  vibrations  of  a  two- 
dimensional  wing  system.  This  study  leads  to  a  more 
comprehensive,  understanding  of  the  (hitter  condition, 
since  it  studies  not  only  the  critical  speed  but  also  the 
approach  to  this  speed.  A  number  of  figures  arc  pre¬ 
sented  that  show  the  nature  of  the  res|»onse  curves  in 
both  one  and  two  degrees  of  freedom.  It  is  perhaps 


worth  mentioning  that  von  Scldippc  (reference  7)  has 
employed  an  experimental  flight  method  for  deter¬ 
mining  tlu1  critical  flutter  speed,  which  is  based  upon 
the  use  of  an  impressed  alternating  exciting  force. 
The  practical  value  of  experiments  of  this  nature  is 
yet  somewhat  doubtful  since  the  flutter  usually  conies 
on  rather  explosively.  In  any  case,  the  theoretical 
results  are  of  interest  because  they  indicate  the  critical 
frequency  as  well  as  the  growth  of  the  maximum 
response  as  the  critical  speed  is  approached. 

METHOD  FOR  ROUTINE  CALCULATION  OF  FLUTTER 
SPEED 

The  cal  ‘illation  of  the  flutter  speed  can  be  reduced 
to  a  routine  procedure  by  the  following  scheme .  Noth¬ 
ing  more  involved  arises  than  the  calculation  of  the 
numerical  value-  of  double  and  triple  determinants. 

(liven  are  a  maximum  of  -even  original  parameter- 
»,  n,  x„.  r.,-,  Sj.  r.  from  which  are  formed  the  I  s  con¬ 
stant*  .1 . \„...  .I,,.  .1,..  etc.,  defined  a-  follows: 
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Those  constants  are  obtained  from  the  original  varia¬ 
bles  and  from  the  T  table  (table  2)  {riven  at  the  end  of 
the  report. 

Another  set  of  quantities  liaa,  etc.,  will  be 
needed ;  their  expressions  are  as  follows. 


1L, 

R.n  + 

/  ..  2f»  71(,  2/ 

2ir  2w  k  2r  f  k:  | 

L, 

/•,.  r 

•Jr  w 

2f«‘  7,..  ln  .,  A 

k  '  2f  2t  -  ) 

R. 

c. 

2<» 

k 

1 ., 

> 

These  si\  quantities  are  derived  from  the  constants 
already  {riven  and  from  two  additional  quantities  /'and 
d.  which  are  functions  of  I  k\  The  quantity  1  /  is. 
in  reality,  the  independent  variable  in  the  problem. 
The  quantities  /'and  (1  occur  in  the  forms  2/’.  2/  it"', 
and  2f»  k:  their  values  are  given  in  table  I  for  different 
values  of  |  k\  In  order  to  facilitate  the  calculation  of 
these  quantities,  the  parts  depending  on  )  l-  are  {riven 
in  tables  .'t  and  I. 

Additional  constants  involved  are  the  frequencies 
51...  !J.».  and  lb.  defined  under  the  scheme  for  each  case, 
and  three  damping  constants '/„.  #/,.  and  (icitcrnlly 
all  these  constants  are  not  simultaneously  needed. 
The  fnui  cases  will  next  be  solved. 

Case  1.  The  problem  is  iriveti  by  two  quadratic 
equations,  for  convenience  referred  to  as  the  “real" 
and  the  “imaginary"  equations.  The  ciH'flicients  of 
each  are  {riven  in  the  calculation  scheme  presenteil  in 
the  follow inir  section.  The  coefficient  of  the  first  term 
in  each  equation  involves  the  constants  </.  and  </,.  the 
coefficients  of  internal  friction  or  structural  dnmpiutr, 
w  hich  are  {riven  as  original  constants.  The  roeflicienl 
of  the  second  term  of  each  equation  involves  the  <7's 
and  the  /f's  and  /' s,  just  defined.  The  constant  term 
in  each  of  the  equations  is  obtained  by  the  schematic 
arrangement  shown  in  the  calculation  scheme:  it  is 


substituted  into  the  equations  and  the  solution,  that 
is,  the  value  of  A’,  is  found.  The  real  equation  usually 
has  two  solutions,  and  the  imaginary  equation  usually 
has  one.  The  values  of  A",  or  preferably  of  ^  .V,  are 
then  plotted  against  ]  k\  and  the  procedure  is  repeated 
until  continuous  curves  representing  the  two  equations 
arc  obtained.  (Attained  judgment  or  the  knowledge 
of  the  solution  of  similar  cases  may  considerably  reduce 
the  lalair  involved  lava  use  it  is  then  possible  to  choose 
reasonable  values  of  I  k  at  the  start.  For  wings  and 
ailerons,  1  /  is  usually  less  than  .1.  very  often  around 
I  or  2.)  The  |mint  of  intersection  of  the  two  curves 
represents  the  flutter  point.  Head  ofr  the  values  of 
A  and  \:k.  The  flutter  speed  is  then  {riven  by  the 
expression 

r,.u;..h  I  1 

'  x*  *x.Y 

Case  2.  The  coefficient  of  the  first  term  in  each  of 
the  two  quadratic  equations  again  involves  the  con¬ 
stants  of  internal  friction  »/.,  ami  ;/»■  The  coefficient  of 
the  second  terms  is  built  up  as  in  case  1.  The  constant 
term  is  built  up  likewise.  Proceed  as  outlined  for 
case  I.  The  critical  speed  i'  then 

u I  I 


where  1  /  and  A  are  the  values  at  the  intersection  point 
of  the  curves  representing  the  real  and  the  imaginary 
equations,  respectively.  There  arc  usually  two  critical 
speeds. 

Case  3.  Case  It  requires  a  more  laborious  calculation 
of  the  constant  terms;  otherwise,  the  procedure  is  the 
same  as  for  cases  1  and  2.  The  flutter  speed  is  {riven  by 

»„«.,&  I  I 

vs  <  v.Y 

Three  degrees  of  freedom.  The  case  of  three  degrees 
of  freedom  requires  the  solution  of  two  third-degree 
equations  in  X.  The  constants  ,,f  the  first ,  the  second, 
and  lilt'  third  terms  are  readily  recognized  as  containing 
only  quantities  already  used  under  eases  1,  2,  and  :t. 
The  expressions  for  the  constant  terms  of  the  two 
equations.  If  and  !>',  involve  three-row  determinants 
but  call  be  obtained  by  strai{rhl forward  calculations 


made  up  from  certain  constants  /#,.  and  //, 
together  with  the  quantities  2/',  2/'  kJ.  and  ’2fl,k.  The 
quantities  //, ,  (\,  ami  llt  are  simple  determinants 

built  up  from  the  constants  etc. 

The  coefficients  of  the  two  equations  must  lie  calcu¬ 
lated  for  a  fixed  value  of  I  k ;  these  coefficients  are  then 


for  each  value  of  1  k.  The  point  or  |>uiuts  of  inter¬ 
section  of  the  two  curves  representing  the  equations 
are  again  representative  of  the  critical  speed,  which  is 
given  by 

r„us„h  I  1 
r  xV  <  v.V 
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CALCULATION  Bl  HEMt 

Case  1  (A,  o).— 

Real  equation: 

Coefficient  of  A"3:  -  '/»'/« 1 

( 'oeffieient  of  A’.‘  ;/«/««•  •  If i/*7'** 

Constant:  .l/|*  . li  *  /fr^.  +  C|^y 

I  iiiittrinurv  ei|imtion:  ' 

Coefficient  of  A’3.'  li»Uj;/.  •  ;/- 1 

Coefficient  of  A’.'  if ■  /«.*  •  if .t//,*;/.  ■  7» 1 

Constant:  ].(/>,  ■ 
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Heal  equation: 

Coefficient  of  A’3:  ifjlf*i  I  ;/.</• 1 

Coefficient  of  A":  ifit/?,*  —  <7t/,* 1  •  Ili'/fu 

—  I  /  —  I  \  Jfi  .’/■ 

Constant:  U*  .1..  /<.  «..*,)  +  C; 

Imatrinnn-  equation: 

Coefficient  of  A‘3:  IJ,if, <</.<  •  ;/. ' 

Coefficient  of  A":  iM/f,-**/.*  •  /,» 1  •  if»(/fMf/»  •  7 m1 

Note  that  when  the  friction  raffle*  nts  #  afr  rt  f**,  a  furl.*  I  *  .^n  »*  rant*  V*1  otil 
«»f  all  term*  In  ait  imaginary  ee|iini*»nv 
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Cate  3  (ii.  A. 

Real  equation: 

(  'ocdiciclll  of  A 
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Imairiniirv  equation 

Coefficient  of  A\ 
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Three  degrees  of  freedom  with  friction. 

Krnl  r«|ii>itiiiii 

<  orflicirnt  of  A Jl.tV.o1 1  ;/„</.  «/<;,  ■ 
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‘rin*  iMjimtimi'  are  urittru  explicitly  fur  1  /.■  1,  that 
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lliv  follow  mt*  <\iiiii|iIi  r«-fi-r»  !•»  ii  tm-lrni  Inter 
iiir|i|mn-  Tin*  parameter*.  wIimIi  weir  (urm-lti’il  l*> 
llir  iiminifiii  liiirr,  mi 


Case  1  i  flexure-torsion 

|{rnl  «-<| mi i  j< >n : 

Tin-  i.if!i.i.iit  nf  A"  i- 
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*1*1  ii*  i-iN'Ilii'ii'iit  i»l  A  i- 
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Tin-  l i’ll  I  ripialioU  l-  t  In'll 


A  mtIhiI  dncrtplion  of  tlir  repmentaliAe  pnniinrlcrs 
used  in  tin  example  I-*  shiv’  den«ity.  «.  itlmiil  2  •*, 
l»>iin,U  |»-r  «|iinrr  foot  per  <- )■•  >r •  I  lenylli  'in  fli. 
«f  iffnn«-ax|s  I* w’nt I* *n  .10  perrenl  "I  tin’  rlinril  fn>m  the 
It-stiltitsT  e>htr.  r»’t»ii’r-«f-tfrin tty  location.  <i»  |>rnrni  «»f 
till'  rln>r»|  In >I||  till’  IrndlliC  eiljre.  aileron  It-nj-l Ii .  olir- 
liflli  of  tlir  total  chord;  halaiii'cil  nili'rnii  >>'i»*ili*r  of 

irrn\ii>  >»f  ntl>‘r»tt  <>t  Inner  axt*  t,  •»’  Tlir  «iriirtural 
datnptne  na’llirirnt*  >/,.  and  »/,  will  lw  Kept  rrm 
>i>rri"i|>»nilme  l»  »  '•airly  fi»  t>>r  It  i«  n»t  ninn-arx 
ti>  -|W!if\  |||,<  chord  Imcth  Jft  and  tin’  t»n*i>miil  frr- 
i|in'tn  a  at.  until  the  final  *trp  The  follow  me  fre- 
«| Ill'll'  A  ratlin.  Il»*w rl  rr .  it r>’  *|ieclfird 


[a-  |  .VJ I  :i  A'  .V70.*i24  It 


Imaeinary  >-«| <i>i t i>>n 

Tin-  riH*  Hi«-»i*nt  i>f  A  i’¬ 


ll!  u  ii  n77n:t  M  ii7am  1  2s4*>4 


Tin*  i'i>ii»tiiiit  term  A//  i- 

:t  :t  ti  7i  11.21  wo  i  i  n.ij'iiii  .iiTss  i  :t.:t ii.*,.*, 7 

Tlir  itnaeimiry  <•<( im t i>»n  i-  tin'll 

1  JMVt.Y  .{,'i(i.V*>7  0 


Tlir  nmt’  <>f  tlir  mil  ri|iuitii>n  nrr  A"  1 . 1 S7  mill 
MM2I.  mul  tlir  root  nf  tlir  imiifin.’irv  eipiation  is 
A'  2..‘»7.t.  ur 


Tlmt  is.  tlir  torsional  frr»|tienry  ii  four  limn  tlir  l*rml- 
irur  frripirnry  mnl  tlir  aileron  fmpirnry  i«  1.22  limn 
tlir  Ih’iiiImut  fmpirnry.  Tlir  r»n«t»nt«  from  wliirli  nrr 
rmii[»wril  nil  tlir  ilrlrrminauti  in  tlir  relriilotinn 
srlirmr  arr  taluilatnl  a<  follows: 


x  A'  l  .OSlt,  4.407,  nod  1.004 

Tlin<’  values  of  ^  A'nrr  plotlivl  apiiosl  1  ,k  in  ftpttrr  1 . 
Tlir  runn  tmnil  by  plotline  tlm  roots  nrr  shown  in  the 
lieurr.  Tlir  intrrsn’tion  is  at  >A  I. *>04.  I  f’  2. 40. 
Tlir  flutter  «|>rril  is  thru 

1  .  2.40 . 
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MECHANISM  OF  FI.ITTER 


In  tho  present  example,  the  chord  2ft  is  12  feet  and  u, 
is  00  (corresponding  to  a  torsional  frequency  of  850 
oyles  per  minute);  bua  is  then  540  feet  per  second  or 
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Fi«.i  kf.  I.  (*um*  t.  Numerical  example.  Ki«;vuk  l.-  t’ase  i.  Numerical 


Tile  i«nds  \  X  of  the  rail  awl  tlu*  imagi¬ 
nary  «  i-*ns  against  l  K. 


example.  Tlu*  mots  \  A" of  the 
real  an<l  tlu*  imaginary  luna¬ 
tions  against  I, it. 


A/,’  -  A/.' 

>1  -Mr1 


«r,  oil  solving  for  A/  - 


about  80S  miles  per  hour.  Hence  the  flutter  speed  is, 
for  this  ease.  507  miles  per  hour.-’ 

Case  2  (flexure  aileron). 

Real  equation: 

The  coefficient  of  A’-’  is 


;j  '0,0012  •  I  0.001S 


The  coefficient  of  A"  is 

0.00 IS (  —  5  •  0.2000)  *•  I  (-0.007028 0.020470) 
0.004808 


•  The  compressibility  correction:  l,»*i  t  l*t*  calculated  lint  tor  s|h*««|  for  the  miom- 
pressihle  flu i<  1  he  r,  ami  lot  tlie  corresponding  s|iee'l  f*»r  the  compressible  fluid  hr  t 
Denote  r,  r  h>  A/,  and  / .  <•  by  M,,  where  <*  is  the  velocity  of  sound.  Then  tsee  Intro- 
duct  inn  i 


— ) 

For  example,  withr,-.V»7m.  p.  h.,  A/*»Sfl7/7Wi-fi.7W,  A/r  =  lMi*«n.  ainl  ),=  lii|in.  p  h. 
Note  that  the  example  given  refers  to  sea  level;  at  altitude,  the  example  should  he 
based  on  another  value  of «  and  an  appropriate  value  of  the  velocity  of  sound. 


The  constants 

A.  0.084001 
.1.  — 0.0S255 

II.  0.00S  1 54 
—0,010510 
C,  =  —  0.024080 
!).  --  — 0.1025S 

The  constant  term  .1  /,"  -  0.07822. 

The  real  equation  is  then 

O.OOlSA'-'-i  0.004S08A'-0. 07822  (I 
Imaginary  equation: 

The  coefficient  of  A’  is 

0.0018X1.078841(0.021 177  -0.0003184)  0.028487 
The  constant  term  M,'  -  —  0.0SS554 
The  imaginary  equation  is  then 

0 .02348 7 .V — 0 .  ( ISN554  0 

The  roots  of  the  real  equation  are  A’  5.182  and 
—  7.85  and  the  root  of  the  imaginary  equation  is  A 
8.778  or  (for  the  positive  roots)  ^  A'  2.270  and  1.044. 
These  values  of  \  A"  are  plotted  against  !  ,/■  in  figure  2. 
The  curves  traced  by  the  roots  are  shown  in  the  figure. 
Since  no  intersection  exists,  this  case  is  stable. 

Case  3  (torsion-aileron). — 

Real  equation: 

'Phe  coefficient  of  A’-  is 


.  8  0.0012 
I  ...>A  ,, 


0.00045 


'  82  0.25  . 

The  eoefiieieut  of  X  is 

l(  0.007028  :  0.020470)  •  0.00045(  1 .285  0.12508) 

0.012807 

The  constants 

._P,  0. 008408  r:i-  -  0.008120 

.  1,  0.021110  (1.00424  1 

lb  0.000700  D,  0.020004 

/f, 0.002000  77,  -0.001474 

Tin*  constant  term  A / /*’  0.021178 

The  real  equation  is  then 

0.0004 5.Y-’  0.0 12807.' V  0.021178  0 

Imaginary  equation: 

The  coefficient  of  A  is 

(0.021177  :  0.0008184)  •  0.00045(0.0  0.07708) 

0.021805 

The  constant  term  MJ  0.080070 

The  imaginary  equation  is  then 

0.02 1805. V  0.080070  0 

The  roots  of  the  real  equation  are  .V  1.507  and  —80.08 
and  the  root  of  the  imaginary  equation  is  X  1.875 
or  (for  the  positive  roots)  \  X  1.252  and  1.178. 
These  values  for  the  \.V  are  plotted  against  1  k  in 
figure  8  'Phe  curve  traced  by  the  roots  is  shown  in  the 
figure .  Since  no  intersection  exists,  this  case  is  also  stable. 

Three  degrees  of  freedom  (flexure-torsion-aileron). — 
Real  equation: 

The  eoefiieieut  of  A’:l  is 

I  •  0,00045  •  0.0001  125 
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The  coefficient  of  A'2  is 

0.00045  ( -  4 .7994 )  +  0 .000 1 1 2 5  ( -  1 .4 1093) 
+  J4  (0.013442) =0.00 1042 


The  coefficient  of  A’  is 

1(— 0.07322) +  0.00045(5. 76524) +  h(— 0.021 173) 


The  constants  are 
R=  -0.032848 
R=  0.077092 
S=  -0.004381 
S= -0.000344 


07592 

7=0.017042 
7=0.028790 
T=0. 103485 
F=0.0 17720 


The  constant  term  /4R=0. 094635 
The  real  equation  is  then 

0.0001 125A'M-0.001042A'2-0.07592A'+0.094635=0 
Imaginary  equation: 

Tlie  coefficient  of  A'2  is 

0.00045(1. 0788) +0.0001 125(0. 82297) +  )<(0.02 1495) 
=  0.005952 

The  coefficient  of  A'  is 

1  (  —  0.0S8554)  +  0.00045  (- 3.30557)  +  ?((-  0.030076 ) 
=  —0.097501 
The  constant  term  Dr=  0.11711 
The  imaginary  equation  is  then 

0.005952A'2— 0.097561  A'+0. 1 17 1  =0 


The  positive  roots  of  the  real  equation  are  A'=  1.270 
and  21.0  and  the  roots  of  the  imaginary  equation  are 
1.302  and  15.08,  or  V-V=  1.120  and  4.58,  and  1.141 
and  3.883.  These  values  of  ■>] X  are  plotted  against 
\',k  in  figure  4.  The  curves  traced  by  the  roots  against 
1/A"  are  shown  in  the  figure.  The  intersection  is  at 
VA'=1.00,  l/A-  =  0.875.  Hence 


;  %  X  2  X  -  j  6o>„  —  0.8266«„ 


For  26=12  feet  and  w„  =  90,  the  flutter  speed  is  304 
miles  per  hour. 

These  examples  have  been  selected  from  several 
listed  under  the  last  part  of  the  billowing  section,  to 
which  the  render  may  refer  for  other  examples,  includ¬ 
ing  the  case  of  an  unbalanced  aileron. 

THEORETICAL  SURVEY  OF  THE  EFFECT  OF  THE 
FLUTTER  PARAMETERS 

The  purpose  of  this  section  is  the  study  of  the  effect 
on  the  critical  speed  of  the  various  independent  varia¬ 
bles.  Although  the  theory  in  itself  permits  the  solution 
of  any  particular  case  without  difficulty,  it  is  somewhat 
difficult  to  obtain  a  perspective  of  the  effects  of  the 
parameters.  Because  of  the  many  variables,  this  survey 
has  been  limited  to  the  magnitudes  and  the  ranges  of 
most  practical  interest.  It  is  realized  that  the  effect  of 
increasing  or  decreasing  a  certain  parameter  is  depend¬ 
ent  on  the  values  chosen  for  the  others.  As  a  mathe¬ 
matical  experiment,  it  is  possible  to  change  one  variable 
ami  to  keep  all  the  others  constant.  With  reference  to 
practical  problems,  however,  the  change  of  one  param¬ 


eter  is  usually  accompanied  by  unavoidable  changes 
in  several  of  the  others.  This  fact  must  be  kept  in 
mind  when  actual  or  proposed  changes  intended  to 
increase  the  flutter  speed  of  airplanes  are  considered. 
This  discussion  is  intended  to  give  only  the  salient 
facts;  the  charts  contain  the  complete  data. 


CASE  I  (FLEXURE-TORSION) 

The  flutter  speed  for  case  1  is  plotted  in  the  coefficient 
form  v,'bua.  In  the  following  graphs,  the  frequency 
ratio  o),,/o)„  is  generally  used  as  abscissa  and  the'critical 
flutter  coefficient  r/bua,  as  ordinate. 

The  graphs  under  each  of  the  following  sections  of 
case  1  are  arranged  in  order  of  decreasing  values  of  k 


Fi'ifKK  3.-  Casr  3.  Niimrrtual  ox- 
ample.  Thu  Tools  >  X  <4  (lie  real 
amt  I  he  imaginary 
against  1A. 


0  /  ^  3 

!tk 


Fk.i  KK  I  Three  «tegrees  of  free, 
•loi.i.  N nine!  teal  example.  The 
roots  ^  X  of  the  real  ami  the 
imaginary  equations  against  \,k. 


starting  with  k  -1/2  (lightest  wings)  and  ending  with 
k  1  2(1  (heaviest).  The  range  of  k  for  present-day 
airplanes  is  approximately  a  1/3  to  k--  1/15.  The 
graphs  are  further  arranged  in  order  of  increasing 
values  of  a,  starting  with  the  smallest  values  of  a 
(stiffness  axis  in  the  most  forward  location).  In  most 
eases,  the  radius  of  gyration  is  kept  at  a  fixed  value 
r.2  1/4. 

Effect  of  center  of  gravity  xa. — The  effect  of  xa  on 
flutter  speed  is  given  in  graph  1  A.  It  may  be  observed 
that  there  is  usually  a  decrease  in  the  critical  speed  ns 
the  frequency  ratio  <■)»/«„  is  increased  from  zero  and 
that  the  curves  tend  to  a  minimum  near  the  frequency 
ratio  wt/u,=  l.  There  arc  cases,  however,  in  which 
the  minimum  critical  speed  lies  at  &>*/&>«=().  The 
transition  takes  place  for  a  certain  small  value  of 


(i  *>t j«f >!*';  /.,=-iu  stable;  p/,/ftw„  *  I.JSh.  <h)  «-  !*»;  a  =  -0.3;  /.,=»  stable;  rn-ta.,-  1.12.  id  *  =  1 1.  n  =  -b.2;  *0.913. 

u|>  $:«;  a=-«.4;.r„=o  stable;  r/,/ftw„  *1.04.  ««  !:<;«  =  -0.3;  n»'nw»*  1.37. 

<f>  «-*s;a«-u.2;rj>/Aw„-|1l2.  *p)  «*  1 •;  a  =  -0.4;  J-«*0  stable;  rvto««2.24. 

(h)  o  *-0.3;  r/j/Aw«*I.SN.  (i)  h;  n  »  -0.2;  rn!bu«~  1.29. 

Graph  I-A  (a  I).— The  effect  of  /*;  the  flutter  coefficient  apainat  the  frequency  ratio;  r*:=  1/4.  Case  1  (A,  o). 
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This  value  is  greater  the  larger  the  values  of  k  (light 
wings).  For  instanee,  when  k=1/4,  a  value  of  /„= 
about  0.1  (graph  I- A  (g))  brings  the  minimum  near  the 


speed  near  unity  frequency  ratio  from  zero  to  infinity. 
As  may  be  observed  later,  structural  damping  will 
greatly  alter  the  shape  of  the  curve  in  this  range. 


origin.  For  «  1/10,  r„  must  be  close  to  zero  (graph 

I-A  (m))  to  cause  transition.  The  transition  is  critical ; 
graph  I-A  (in)  shows  that  a  2..!>-percent  change  in  the 
position  of  the  center  of  gravity  chnnges  the  flutter 


The  range  of  most  practical  importance  is,  however, 
the  neighborhood  of  the  zero  frequency  ratio.  (lor 
wings,  the  ratio  is  approximately  1/4.)  In  this  range, 
the  parameter  of  greatest  significance  is  really  the  com- 


MECHANISM  OF  FLUTTER 


15 


bination  ^-\-a-\-xa.  1»  other  words,  the  flutter  speed  light  wings  (k  I  d  to  I  d),  (iraph  I  ('  hi;  gives  very 

is  very  nearly  a  function  of  the  locution  of  the  center  of  normal  values  of  the  parameters  as  used  for  most  wings, 
gravity  with  respect  to  the  forward  quarter-chord  posi-  These*  curves,  for  a  tjm.n  wing,  may  he  taken  to  give  the 
tion  and  not  of  the  distance  relative  to  the  stillness  effect  of  altitude.  Note  that,  for  a  given  wing  with 
axis.  Graph  I-B  (a)  shows  clearly  that  the  value  of  a  k=1/5  at  sen  level,  k  becomes  1/10  at  approximately 


(iraph  I- 1)  i.a  f'.-  The  tiltri  uf  Mniciiir:il  friction;  !).«•  Iliiltci  nu  Hicii-rii  iii>:ttn^(  ihc  ■!!«•>  rain*,  a  |jo:«  -n .2:i.  n.l.  Cast*  l 

actually  has  no  inlluenee  on  the  (lutter  speed.  Outside  j  15,000  feet,  with  a  resulting  increase  in  the  flutter 
of  this  range,  that  is,  for  larger  values  of  the  rela-  speed  under  normal  circumstances.  For  the  case  with 

tionship  is  less  simple,  (iraph  1-11  (b)  shows  the  s„-  0.2  given  in  graph  1  C  pi),  the  increase  in  tin 

dependency  on  the  center-of-gravitv  location  lor  various  Mutter  coefficient  is  from  1 .0  to  1 .05,  or  about  20  percent, 

positions  of  the  stiffness  axis  a.  For  a  constant  xa,  j  It  is  possible  that,  for  very  light  wings,  the  flutter 

that  is,  for  a  constant  distance  between  the  stiffness  I  speed  might  decrease  with  altitude  until  a  certain 


•:»  *  I*'  »  '  U  *  i 

(irupli  I  K  it  rt.  Tin*  rlTcti  «»1  rixiiii.*  *'l  p\r;ito<ii.  I  h«-  Ihilict  iurllirtriil  ti^uin.sl  thr  frfjiH'mj  rat  i<  <;  n  - 1'  j..  n.l.  Cum*  1  i  A.  <*). 

axis  and  the  center  of  gravity,  the  flutter  speed  is  |  altitude  i<  readied.  For  high  values  of  I  k  (heavy 
increased  us  the  stiffness  axis  (and  center  of  gravity)  i  wings),  the  flutter  speed  increases  nearly  ns  the  square 
is  moved  forward.  ;  root  of  the  wing  density,  1  *. 

(iraph  I  ('  shows  the  flutter  coefficient  plotted  against  I  Effect  of  structural  friction  </„,  Graph  1  I)  is 
I  k.  The  normal  range  of  wings  is  included  in  the  1  intended  to  show  the  clfeet  of  the  structural  Irietion  on 
diagram  (the  heaviest  wings  to  the  right).  The  ilia-  the  critical  speed.  As  the  eoeflicients  of  friction  are 
grains  arc  arranged  in  order  of  increasing  values  of  increased,  there  is  a  definite  tendency  for  the  often 
(«*'&.'„)'  and  of  rt.  An  interesting  result  is  the  existence  pronounced  minimum  flutter  speed  near  on,  u!„~l.(l  to 
of  a  minimum  critical  speed  that  falls  in  the  range  of  disappear  and  to  produce  response  curves  of  the  type 


h* 
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obtained  for  negative  value  of  sa-  In  the  range  of  most  (Values  in  the  preceding  table  are  given  in  relation  to 
practical  interest  (u/iiaSO),  the  torsional  friction  is  the  value  for  1/5,  r«2 ----- 1  which  is  the  case  of 
the  more  important.  lowest  wing  density  and  smallest  radius  of  gyration.) 

Effect  of  radius  of  gyration  ra. — Graph  I-K  is  nr-  The  speed  corresponding  to  given  stiffness  drops  if  any 
ranged  in  conventional  order.  Note  that  the  flutter  mass  is  added  so  that  /«,  the  density  I  k,  or  both  are 
coefficient  in  the  low  a ihlua  range  increases  with  in-  increased.  Hence,  any  mass  added  not  for  the  purpose 
crease  in  the  radius  of  gyration.  This  increase  in  the  of  increasing  the  stiffness  or  moving  the  center  of 
flutter  coefficient  does  not  necessarily  correspond  to  an  gravity  forward  is  detrimental. 

increase  in  the  flutter  speed:  it  <lms  if  the  torsional  Flutter  frequency.— The  flutter  frequency  is  shown  in 

frequency  is  kept  constant.  If  the  stiffiirxs  is  graph  1  K.  It  is  seen,  for  instance,  that  for  small 


fm  *  •  a  -0.4.  <!•>  —  u.2. 

«  >  m:  a  -  -0.4.  t«S>  k  1 1»,  a  -  —  0.2. 

Otaph  1  Kg*  *1 FlmiiT  frr«|Ufnt>  rutin  tut  <h-|Manilt'iif  on/.,  auain^t  frrqiU'nry Taliojr.,*'-  )t.  rust*  1 

kept  constant,  which  means  that  w„  is  decreased  as  values  of  wh  u'a,  the  flutter  frequency  is  around  titl 

1  / rm ,  flic  flutter  speed  is  actually  decreased,  as  is  shown  percent  of  the  torsional  frequency  co„;  for  higher  values 

in  the  following  table.  of  the  flexural  frequency,  the  flutter  frequency  ap- 

FI.ITTKR  SI’FFli  FOR  CONSTANT  TORSIONAL  proaches  or  exceeds  the  torsional.  This  graph  is 

STIFFNESS  primarily  of  interest  in  connection  with  experimental 

l«*=-(i.2,/.-o.i,  («,/«.)>- «i  flut  ter  research. 

Coupling  factor  {. — Consider  a  two-dimensional  case 
of  flutter  in  which  only  a  part  of  the  total  length  of  the 
(infinitely  long)  wing  is  given  the  second  degree  of 
freedom.  This  arrangement,  because  of  the  deficient 
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N 

_ 

1,5 
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1/5  ! 
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07. « 
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coupling,  exhibits  a  higher  critical  speed.  Call  the. 
fraction  having  both  degrees  of  freedom,  f.  The  rcsulis 
are  shown  for  several  values  of  (  in  graph  1  (i. 

Divergence  velocity  and  approximate  flutter  for¬ 
mula. — It  can  he  shown  that  the  divergence  velocity 
may  be  expressed  in  nondiinensional  form  as 


to  oil  percent  of  the  total  chord  represents  a  stabilizer- 
elevator  or  a  fin-rudder  combination.  Several  values 
of  r,i  and  //  and  of  the  damping  coefficients  </,,  and 
have  been  included. 

It  should  he  mentioned  that  ordinarily,  as  shown  in 
reference  1.  case  2  differs  hesicelh  from  case  1  bv  the 


The  divergence  velocity  r,,  hu„  is  given  in  graphs 
I  A.  This  velocity  is  usually  higher  than  the  flutter 
velocity. 

An  empirical  expression,  which  is  useful  in  quickly 
obtaining  the  order  of  magnitude  of  the  flutter  speed 
for  small  values  of  oia  and  which  appears  to  hold  very 
well  for  heavv  wings  (with  k  1  11))  is  given  bv 


Graph  I  B  fa)  shows  the  curve  obtained  from  the 
empirical  expression  (dashed)  and  a  curve  based  on  the 
exact  values  (in  full  lines). 

CASK  a  'FLEX I  KK-Alt.KKON 

The  flutter  coefficient  for  case  2  is  r  Im , .  'Flic  fre¬ 
quency  ratio  a).,  is  ordinarily  used  as  abscissa.  The 
graphs  are  again  arranged  in  order  of  increasing  wing 
density.  Two  values  of  the  location  of  the  aileron 
hinge  axis  r  have  been  included.  The  first  value,  r 
or  the  aileron  chord  equal  to  2.1  percent  of  the  total 
chord,  is  intended  to  represent  a  wing-aileron  combina¬ 
tion;  the  second  value.  r~  0.  or  the  aileron  chord  equal 


r  f  8  ‘0  >C  <3 

(1.  lit.  .  fT.  .1  if  flu  >  ■  Uflit..-  *  i>  «•  f  lit.  Muri.r.  •(*.■! 
ir.  1 1 1 1« j h  >  r»t  i".  s  'i  “  j  1  "  : 

existence  of  a  (Inthr  ramji  extending  between  a  lower 
and  an  upper  flutter  speed.  This  range  of  flutter  can 
be  reduced  or  eliminated  by  various  means,  h  i-  im¬ 
portant  also  to  notice  that,  beyond  a  certain  value  of 
the  frequency  ratio  w. ,  a >,,,  in  fact,  fora  value  slight  I  v 
greater  than  unity,  no  critical  speed  cxi-ts.  since  the 
critical  area  docs  not  evi.-nd  much  beyond  tbis  point, 
i  The  reduction  of  the  c<  -of-gravity  distance  from  the 


J~j  -  /JO 


9,  -9-  0, 


\r.‘  /SC 


Coup  og  fcctor.  i'S  | 
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hinge  lias  t lit*  effect  of  reducing  and  finally  eliminating 
the  critical  flutter  area.  Internal  dumping  shows  the 
same  general  effect.  The  fact  that  the  aileron  extends 
effectively  over  a  shorter  length  is  theoretically  ex¬ 
pressed  by  a  “coupling  factor”  £.  which  is  the  length 
of  the  aileron  divided  by  the  total  length  of  the  wing 
executing  deflection.  The  effect  of  £  is  shown  in  some 
of  the  graphs. 

Effect  of  frequency  u>„  (r  *4). — - Graphs  )|  A  tat 
and  (b)  show  the  effect  of  varying  x.f  in  reducing  the 
critical  area.  The  effect  of  damping  is  shown  in  graph 
11  A  (c)and.  finally,  the  effect  of  r/  in  graph  11  A  (di. 

Effect  of  center  of  gravity  x.,  ic  Graph  II  M 

shows  the  flutter  coefficient  against  the  eenter-of- 
gravity  distance  x.,.  giving,  for  two  values  of  *.  the 


1  lie  defined  for  each  \alue  of  the  frequency  ratio.  Il  is 
necessary  then  to  choose  the  largest  frequency  ratio  or 
the  smallest  unbalance,  then  to  calculate  the  other 
value,  and  finally  to  choose  the  most  practical  combi¬ 
nation,  using  a  margin  of  safety. 

Effect  of  radius  of  gyration  ra  n  Graph  II  -(' 

|  shows,  for  a  typical  wing-aileron  case,  the  effect  of 
changing  the  radius  of  gyration  for  various  values  of 
the  frequency  ratio. 

Effect  of  frequency  w., ir  (ii.  In  tin1  preceding 
graphs,  tlie  hinge  axis  was  at  <  '...  Graphs  II  l>. 

II  K,  and  II  V  show  the  results  for  c  tl.  The  cuncs 
j  are  arranged  in  order  and  show  the  effect  of  /■„,  ,,  . 
!  and  ;ij  for  k  and  'in.  One  curve  is  also  included 
.  for  k  (graph  1 1  1)  (d  l) . 


ort, i>i,  u  a  .,  <n  Klim,  i  tin  llirirtif  /  -  f..r  \  tri'-tjs  in  <;ii.  tn  v  r.t*ii-«.  I’  ts*  j l.fr 


effect  of  varying  the  frequency  ratio  u,  at  three 
values  of  r j-.  Note  that  for  large  (beyond  normal 
range)  the  type  of  flutter  reverts  to  that  of  case  1 ;  that 
is,  the  upper  flutter  speed  becomes  infinite  for  a  certain 
value  of  /> 

It  is  important  to  notice,  by  considering  each  curve 
in  this  figure,  that  must  be  decreased  below  a  certain 
value,  which  is  rather  critical,  in  order  to  avoid  flutter. 
If  is  larger  than  this  value,  the  lower  flutter  speed 
remains  at  a  virtually  constant,  small  value.  The 
frequency  ratio  exhibits  a  similar  effect;  that  is,  flutter 
is  eliminated  beyond  a  certain  frequency  ratio  often 
greater  than  unity,  whereas  for  smaller  ratios,  the 
lower  flutter  speed  remains  at  a  low,  nearly  constant 
value.  In  other  words,  a  critical  frequency  ratio  can 
be  defined  for  each  value  of  tin-  unbalance  and,  in¬ 
versely,  a  definite  critical  value  of  the  unbalance  can 


Effect  of  center  of  gravity  j\i  ic  ()>.  The  figure-, 
an-  given  in  graph  II  K.  arranged  as  usual. 
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Effect  of  coupling  factor  (c  ()i  in  graph  ll-K  the 
effect  of  the  coupling  factor  £  is  shown  for  an  extreme 
case  of  unbalance  fx*  large).  The  superimposed  effect 
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of  damping  i-.  'Iiowii  for  tin-  /mi  (roiui'iioy  ratio. 
Notice  liou  llir  t'ouplinp  ftirlor  gradually 

eliminate*  llir  Hotter  u mi. 

(AMK  1  T0MSI0S.»II.K«0\ 

Three  graphs.  Ill  A.  Ill  O.iiikIHI  t  *.  are  pn-srutnl. 
Tlirri'  is  ti  similarity  lo  nisi-  2.  <ini|ili  111  A  show* 
Iiow  tin1  internal  damping  merm'i*-  the  lower  Muller 

‘  "*  '  J  / 


'  hnlitnri'  I/,-  O'  i«  1 1 << -  — ;i i\  |al  iluiiuiule  It 1 1 1 1  •  . 

!  hor  light  win”',  ili”  iir.i-i  i'  It"  | irm ion i ii*  ■,  I  ami 
i,  0  i«  ii'imlly  'iiHii  ii'iil  1 1  tuny  Im  observed  from 
|  (In-  original  -el  of  «*«| mi t  loti'  llial  I r  in*  hahilii'e  agmii'l 
rot  ii  I  ion  iniplic'  /  j-  ir  a  /,  o. 

Timrr:  iiri.m»  or  nttriHiM 

In  order  lo  familial  i/.e  (lie  reader  wiili  the  <  oiupleic 
ni'C  of  din*)'  ih*gree»  of  fri'i'ilmii  iiml  il'  rehilioti'hip  lo 
|  llir  ill  tee  'iiIh-ii')-'.  it  ')•!  of  ivpieal  1  io-ii »«*'  i'  'In iw  1 1 
Tin-  I'oii'imii *  iimsI  mi-  dn'  *iiini'  n»  ill . in  dn'  iiiiim  ii- 

ral  l '  V  it  1 1 1 1  >l<-  (|1.  SI  willl  'OHIO  llililil  loll'  (  ’ll')' 


•  ll(  till  f  h  ffn  »  •  f  i  \»*  »  t  l  |V» 


' jit'i'i I .  (iitipli  III  0  ri'|»rt"*riii'  iluiu  lakiii  (min  mi 
ni-li:iil  i'll'”  of  n  I i^li t  win”  with  n  'iiuillor  nili  ion. 
\olr  llir  sinking  similarity  to  i'll')'  2  For  dir  vain, 
/.i  (i  imi'.iV  lenmplelely  uniniliiiii'i'tl  aileron', 
in  ii  't  In-  greater  llimi  on  to  avoid  llutii-r:  for  tin-  mo;, 
nortiiiil  Milnr  s, t  IMHI2,  u'.,  ui„  nr, si  only  hr  >o  |. 

Tlir  lluttrr  iirrn  i«  ,Tuiiiiiiilr,l  liy  ir, hieing  j,  to  a 
'lightly  smaller  v nlur. 

Case  :t  itor>ioii-nilrroid  i'  probably  of  I,-"  pimticnl 
illl|iorttl»rr  lirrmi'r  llir  rlilllilintioli  of  llllttrr  for  ril'f 


«.,.|  '  III  1*  I  •  • 

•T  *  '  .1 

I  i'  'lioun  in  ligure  I  mnlrl  llir  noiii<-n<al  r\aln|i|r 
Till-  flilllrr  ••tM'lIirirlil  f#to'„  I  ''I-' 

I  "a*r  2  I*  'Iiowii  m  lignir  ", ;  ,  arli  part  of  dir  figure 
refers  lo  illlTrrrlll  roinlima  I  loll'  ol  .  ali«l  ^  ,  N  o 
ft  II 1 1  cr  ornii'  for  dll'  ,  olnliilialioli'  -Iiowii  III  tiglll”' 
*i  a  ami  .*>  It  Ims ai M ',*  of  llir  lialamnl  ailrroti  ami 
liolir  hi  tigiilv  o  <1  Ins  ail','  of  llir  large  all, 'ton  frr- 
iplrlirv  Foi  llir  roinl, million  shown  m  ligut'r  o  ‘r'„ 
tlirir  >  a  lioimal  tahgr  of  flultrr  with  two  finder  |H,iti1s 
'how  ii 

< i'  'Iiowii  ill  figure  « » ;  rarh  pari  refer-.  ir'prc- 
tively,  to  tlir  Milnr  ailrion  parameters  used  m  rii'r  If, 
i  Note  that  ,»• ,  I-  in  al!  •  iws  Tlir  roinl  ana  I  ion* 
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« »r  «|»t*  III  \  PfT*ft  iiMrufh.n  it*  ID,  »  pr,  f  f.  flii*  l»  f  « > «  fT,t  «  rd  vv'sf  'r* 
•Hi*  nr>  Mif>;  I.  N-  1  %.tf  •  .  *  *  ;,«•  f  *«.  : 

2  ordinarily  rsrludr*  tlir  possibility  of  lluttrr  in  rase 
It;  hut  it  is  notr, I  tlint.  in  onlrr  to  rliniinnlr  ninss 
coupling  in  thr  tnrsion-Hilcnm  rns,*.  a  roiuph-tr  balance 
of  llir  nilrron  iii  tlir  ordinary  srnsr  ir,  tli  i»  not  ipiitr 
siiflirimt.  It  is  actually  found  in  thr  ,',«*«•  of  n  liravv 
wing  and  no  intrrnnl  friction  (with  rt  0,  that  tlir 
lluttrr  speed  is  low,  particularly  near  w.,  —  u.  Kvrn 
n  slight  amount  of  friction,  however,  is  suHicirul  to 
cancel  the  cause  of  this  flutter,  (irapli  III  <’  (fairly 
heavy  winp)  shows  that,  for  no  friction,  a  small  over- 
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shown  in  liirurc'  rt  m »  and  ii  ih>  nn>  a>:iiin  stahlc  Urcausc 
of  the  ailcnui  mass  hnlaucc.  For  the  arrangement 
shown  in  liirurc  ti  fd>.  the  ailcnui  fnspicncy  is  not  hipli 
cnoiiph  to  pn'vcut  flutter  as  it  did  in  case  2.  (Condi¬ 
tions  an-  s|j||  worse  for  the  combination  shown  in 
figun-  ii  («•  i. 
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Fur  tin*  him*  «»f  ilirii'  ili*t*rt-fi  uf  fri*-i,iti*iii .  Ii^iiii*  i 
shows  lln-  n—ull-.  arranged  n  i  I  In*  -mm*  order  u*  under 
i'in*i  ’  ami  !  K<>r  tin*  riiiiiliiimii  (.'inn  in  li”iin*» 
7  mi  ninl  7  tin.  Holler  existed  only  in  ease  1.  The 
Hotter  |*oint  shown  is  therefore  tnllv  ran*  I 

lln IIit.  Thi*  value  nf  tin*  Holler  riM*llirii*nl .  however. 


TRANSITION'  TO  I'll  KKK-IHMKN  SION  At.  FU'TTKR 
PROBLEMS 

Till-  previous  theory  It-Iji  I  <**.  In  tivu-dilliensinliul 
IIiiIIit  ami.  'irii  ilv,  in  a  vv in"  nf  infinite  lentil  Ii.  The 
second  IC'I  riel  inn  i'  iml  Very  lluuhlc'onie.  the  aspect  - 
t  n  I  in.  nr  -pan.  elferl  lieinn  relatively  lit  I  iliijii  irl  nut  ami 
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i  »  i  ...  -  >  ■ .  -  .  f.  . 

IliH  ailoallv  fnnn  ll~  ea«e  I  value  nf  I  'il'J  In 

■  >.7n  ami  I*  v.’.'i,  respectively. 

hnr  t l*i*  arraniirment  shown  m  lignin  7  ul  .  Holler 
exists  in  eases  I  mill  •<  Here  I  lie  range'  completely 
merge.  indicating  'lability  al  only  mi  Inn  »|ieeil 
Holler  exist'  ill  all  three  en»e«  fur  the  eniiiliiiialinii 
shown  in  ligure  7  ie  The  ea-e  J  Muller  ran  lie  rec¬ 
ognized.  a  I  mil' I  unchanged.  while  again  the  Ilullei 
ranges  nf  m«r  |  ami  i*u n*  :i  have  merged.  a<  in  figure 
7  fill. 

Figure  s  has  lieen  iiiehuleil  In  'Imw  dial  then*  is  ii 
enii'iilerahle  lowering  of  Hie  Hotter  «|iei*il  fur  low  value' 
of  the  aileron  frequency  even  though  the  aileron  is 
halnneeil.  Thi'  enmlitiou  i'  prohahly  not  nf  primary 
concern  l<ero'*'<*  n  small  amount  nf  frieliun.  particularly 
lf#,  w  ill  re'lnr,  , lie  flutter  s|H*ed  to  its  full  Iran*  1 1  value 
ll  i'  In  lie  imlefl.  however.  I  lull  a  slight  nverhalanee 
(/,  O'  may  la*  desirable. 


I  »  .  I  .1  - . 

■  '  ■  . . . . 1 .1.  i  v..  ,  Il'.l.  fie  ). 

by  nn  lueaii'  a~  great  a*  the  U'poet -la  I  in  eMeel  iismi- 
eialeil  with  stationary  flows  ll  max  he  disregarded 
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ami  Ineilly  eousiilereil  as  a  safely  faelor.  since  nn  air 
'|M*e<l  of  ihr  order  of  a  few  percent  more  Ilian  that  in 
two-dimensional  Mow  is  necessary  lo  cause  flutler. 
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C'onsider  the  ease  of  a  rectangular  cantilever  wing. 
Some  authors  have  attempted  a  solution  on  the  assump¬ 
tion  that  the  response  curves  in  torsion  and  deflection 
under  normal  conditions  (zero  air  speed)  may  he  used 
in  the  flutter  theory.  It  is  contended  that  this  assump¬ 
tion  is  false.  Several  rather  interestin';  experimental 
results  will  he  presented  in  the  next  section,  which 
show  directly  and  indirectly  that  the  modes  in  flutter 
differ  radically  from  the  ordinary  ones.  The  followin'; 
<|iicstions  arise:  (1)  Does  the  win*;  flutter  in  the  first, 
second,  or  third,  etc..  bending  “mode"?  (2)  Are  these 
modes  in  any  way  related  to  the  ordinary  types  of 
vibration  modes? 

Consider  first  the  case  of  a  very  high  bending  mode. 
It  is  useful  to  consider  the  two-dimensional  case  as 
representing  the  "averages’’  of  parameters  and  variables 
of  thi>  three-dimensional  case.  The  variable  h  now 
appearing  in  the  (two-dimensional)  equations  refers 
to  the  “average"  li  which  approaches  zero  even  though 
the  local  h  in  the  loops  is  very  large.  It  is,  furthermore, 
evident  that  the  average  curvatures  are  greater,  the 
greater  the  mode.  Both  these  conditions  are  equiva¬ 
lent  to  a  large  coefficient  of  internal  damping,  since 
the  work  lost  per  cycle  refers  to  a  very  small  average  h. 

v  y  y  y 


Kb.'  uk  Si-hwiitiiic  lit* i ire  f««r  hiiiher*«.nler  hemline  m* •«!«•  of  cantilever  vim: 
restraint  by  wires  arvl  with  deflect  it  nis  in  phase. 

The  coefficient,  in  fact,  rapidly  approaches  an  infinite 
value  as  the  number  n  of  the  mode  is  increased.  It  is 
probable  that  the  second  “flutter”  mode  involves  a 
coefficient  of  damping  20  to  off  times  larger  than  that 
of  the  first. 

A  study  of  the  graphs  with  specific  reference  to  the 
effect  of  damping  shows  that  higher  flutter  modes  can 
be  expected  only  under  very  special  circumstances. 
This  fact  does  not  mean  that  flutter  occurs  in  the  low¬ 
est  (zero  airspeed)  bending  mode. 

The  bending  frequency  in  flutter  of  a  cantilever  beam 
is  determined  by  a  certain  minimum  condition.  The 
wing  will,  of  course,  flutter  at  the  lowest  speed  possible. 
It  will,  therefore,  not  assume  its  lowest  (stationary) 
bending  mode  but  will  tend  to  assume  a  mode  of  a 
higher  frequency.  Since  this  higher  frequency  tends 
to  uncouple  the  h  degree  of  freedom,  the  actual  response 
ordinarily  happens  to  be  a  cross  between  the  first  and 
the  second  modes.  Large  internal  friction  will  tend  to 
push  the  response  closer  to  the  first  mode.  The  result 
is  u  flutter  speed  distinctly  lower  than  that  calculated 
on  the  basis  of  the  frequency  of  the  ordinary  funda¬ 
mental  bending  mode.  The  flutter  speed  calculated 
by  using  the  lowest  bending  frequency  is  too  favorable. 
In  the  case  of  wings  of  small  internal  friction  (solid 
metal  wings),  the  actual  flutter  speed  is  only  about  0.0, 
the  speed  calculated  on  the  basis  of  the  lowest  bending 


frequency.  In  the  case  of  conventional  wings,  the 
error  is  apparently  in  the  order  of  only  1  or  2  percent,  a 
fortunate  coincidence  because  it  permits  the  use  of  a 
small  experimental-empirical  correction.  This  point 
of  view  is  in  harmony  with  the  Rayleigh  principle' 
which  states  that  any  response  function  whatsoever 
corresponds  to  a  frequency  higher  than  that  of  the 
fundamental. 

To  recapitulate:  Tin  bnnliny  fritpnncij  inm/rel  in 
thefluttir  of  a  cantihnr  winy  is  ijriatir  than  that  obsinnl 
at  zero  air  spiel;  tin  in  on  so,  tin  hnnr  tin  inti  null 
(lam /liny  of  the  winy  struct  a  re. 

This  interesting  phenomenon  is  demonstrated  by  the 
photograph  of  the  flutter  of  a  uniform  cantilever  beam 
(see  fig.  IS)  presented  in  the  experimental  section, 
which  shows  that  the  maximum  amplitude  is  not  at 
the  tip  but  is  rather  close  to  the  middle  of  the  span. 

Another  very  convincing  experimental  proof  of  this 
phenomenon,  given  in  the  experimental  section  of  this 
paper,  is  that  a  counterweight  at  the  tip  section  in 
front  of  the  center  of  gravity  actually  Imres  the 
flutter  speed  of  a  uniform  cantilever  wing.  For  a 
relatively  small  counterweight,  the  tip  section  is  beyond 
a  node  in  the  h  curve.  In  this  same  connection, 
another  rather  remarkable  experiment  was  made:  A 
cantilever  wing  flutters  at  about  200  miles  per  hour. 
The  point  where  the  node  of  the  second  bending  mode 
(at  zero  air  speed)  intersects  the  torsional  stiffness 
axis  was  fixed  by  connecting  this  point  by  wires  to  the 
tunnel  walls.  The  wing  subsequently  fluttered  at 
!.">()  miles  per  hour.  The  flutter  stopped  when  the  wire 
broke!  The  explanation  is  that  the  bracing  wires 
“couple”  a  bending  mode  that  was  previously  entirely 
“uncoupled."  It  should  be  noted  that  the  frequency 
actually  involved  in  this  flutter  is  again  in  excess 
of  that  of  the  second  bending  mode  (at  zero  air  speed  i : 
large  forces  are  therefore  transmitted  through  the  wire 
supports  from  the  walls. 

In  order  to  illustrate  more  convincingly  how  the  sup¬ 
port  wires  lower  the  flutter  speed,  reference  is  made  to 
figure  !>,  which  shows  a  high-order  bending  mode  of  a 
wing.  If  this  bending  frequency  is  about  equal  to  the 
torsional  frequency,  the  lowest  flutter  speed  is  obtained. 
When  the  support  wires  are  removed,  the  wing  will  tend 
to  vibrate  about  a  fixed  mass  center  line,  with  the  result 
that  the  average  h  deflection  becomes  zero  and  all  h 
couplings  disappear.  The  a  moments  and  the  It  forces 
transmitted  to  the  support  are  good  measures  of  what 
may  be  called  the  effective  values  of  «  and  h  when  the 
two-dimensional  theory  is  applied  to  three-dimensional 
cases.  For  instance,  the  transmittal  of  a  small  It  force 
to  the  support  indicates  that  the  positive  and  the  nega¬ 
tive  li  values  very  nearly  cancel.  The  h  effect,  although 
locally  large,  m«y  very  nearly  cancel  itself.  This  fact 
does  not  prevent  the  use  of  a  certain  (small)  average  or 
effective  b  in  the  calculations.  With  no  internal  fric¬ 
tion,  the  flutter  speed  is  not  changed.  As  was  pointed 
out  before,  the  use  of  the  small  effective  h  for  higher 
modes  is,  in  reality,  equivalent  to  employing  a  greatly 
increased  coefficient  of  internal  friction. 
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Tliis  discussion  and  associated  experiments  lead  to  the 
important  conclusion:  A  bracing  wire  may  lower  the 
critical  speed  of  a  cantilever  wing  or  fin.  It  usually  does 
lower  the  critical  sliced  when  the  internal  damping  in 
the  structure  is  low. 

For  a  rectangular  cantilever  wing,  there  is  no  difficulty 
In  regard  to  the  other  parameters.  "With  the  bending 
modes  considered  known,  the  variables  a ,  0,  and  It  were 
given  simply  as  average  values  and  used  in  the  two- 
dimensional  solution.  Inspection  shows  that  the  flutter 
speed  of  a  uniform  cantilever  wing  is  essentially  that  of 
the  two-dimensional  case  involving  the  same  parameters 
and  the  proper  frequency  ratio. 

A  cantilever  wing  of  normal  tapered  shape  will  next 
be  considered.  It  is  assumed  that  there  is  a  similarity 
in  construction  along  the  span  for  each  cross  section. 
The  mass  is  put  equal  to  a  constant  times  the  square  of 
the  chord:  static  moments,  to  the  third  power:  and  the 
moments  of  inertia,  to  the  fourth  power  of  the  chord. 
Further,  the  air  force  is  proportional  to  the  chord  and 
the  acting  moments  are  proportional  to  the  second 
power. 

Various  weight  factors  of  the  form  (b  hr)’‘J(j)il.r  arc- 
obtained,  where  J\.r)  is  a  weighted  wing  parameter  and 
./■  is  measured  along  the  span.  If  the  reference  section 
is  chosen  in  such  a  location  that,  for  a  particular  //  — - m, 

=  JVw* 

which  is  always  possible,  then 

In  other  words,  the  proper  choice  of  a  reference  section 
renders  the  weight  factors  of  approximately  equal  magni¬ 
tudes.  If  the  reference  section  is  taken  too  close  to  the 
tip  of  the  wing,  there  will  lie  a  certain  positive  correc¬ 
tion;  if  chosen  too  close  inboard,  there  will  be  a  nega¬ 
tive  correction.  The  correct  value  is  thus  virtually  con¬ 
fined  between  definite  limits.  The  most  representative 
section  will  lie  close  to  the  three-quarter  semispan 
location. 

In  the  two-dimensional  case,  the  length  along  the  span 
is  considered  to  be  equal  to  unity  and  this  unity  is 
treated  as  being  large  as  far  as  span  effects  are  con¬ 
cerned.  If  the  length  is  different  for  the  two  variables 
considered,  a  slight  modification  of  the  theory  is  neces¬ 
sary.  Knelt  length  is  considered  to  be  long  enough  to 
permit  disregard  of  aspect-ratio  corrections  for  the  air 
forces. 

This  sort  of  consideration  is  of  interest  chiefly  in  the 
case  of  ailerons  and  tail  surfaces.  The  equation  giving 
the  equilibrium  of  the  ailerons  refers  only  to  the  length 
of  the  aileron.  The  included  area  of  the  It  curve  is 
sometimes  a  small  fraction  of  the  total  area  under  the 
h  curve.  This  fraction  will  be  called  £. 


The  solution  of  the  deflection-aileron  case  is  given  in 
reference  1  by 

.  j  _  H  bj  ' ~  1 1 r>3  H  bh~ r  hh  ,, 

Iic3  '-  ilc9  Rc,A  iIrh 

and  with  the  effect  of  £: 

m  1 _ Hti3  1 1  b,1  £  (  II  btl  ~i~  I  I  hh)  ,. 

"" "\Rn>+  U'slU-Hr* 

It  is  noticed  that  the  factor  £<(1  describes  a  certain 
uncoupling  of  the  system.  The  calculation  of  flutter 
speed  can  be  performed  for  any  coupling  factor  £. 
Again  it  should  bo  remembered  that  the  free-vibration 
modes  arc  not  identical  with  the  flutter  modes.  A 
tendency  exists  for  £  to  approach  unity  since  the 
aileron  forces  the  motion  of  the  wing. 

EXPERIMENTAL  FLUTTER  RESEARCH 

GENERAL 

The  purpose  of  the  experimental  research  was,  first, 
to  check  the  theory  as  regards  accuracy  and.  second,  to 
provide  a  basis  for  an  understanding  of  problems  met 
with  in  .airplanes. 

These  tests,  about  one  hundred  in  all.  were  conducted 
in  the  X.  A.  A.  8-foot  high-speed  tunnel.  (See 
fig.  10.)  In  order  to  protect  the  propeller,  a  heavy 
wire  screen  was  inserted  in  the  test  section  immediately 
behind  the  flutter  model.  For  convenience,  models 
having  a  flutter  speed  below  :t00  miles  per  hour  were 
tested. 

The  procedure  followed  was  to  increase  the  tunnel 
speed  slowly  until  flutter  appeared.  If  the  flutter  was 
of  a  violent  type,  the  load  was  immediately  dropped  to 
save  the  model.  In  the  tests  on  ailerons,  the  lower 
branch  of  the  flutter  curve  was  similarly  obtained.  The 
upper  end  of  the  range  was  obtained  by  the  following 
method:  The  aileron  was  kept  in  place  by  restraining 
wires  attached  to  its  rear  end  and  running  across  the 
tunnel.  By  manual  operation  of  the  wires  from  flu1 
outside,  the  arrangement  could  he  conveyed  through 
the.  dangerous  range;  on  slackening  the  wires,  the 
operator  would  receive  indication  of  incipient  flutter 
until  the  speed  had  increased  above  the  dangerous 
range.  When  the  upper  stable  region  had  boon 
reached,  the  wires  were  completely  released  and  the 
conventional  flutter-test  procedure  was  reversed:  that 
is,  the  tunnel  speed  was  slowly  decreased  until  the 
violent  flutter  appeared.  The  restraining  wires  were 
then  immediately  tightened,  and  the  speed  was  noted. 
The  effect  of  the  very  fine  wire  was  shown  to  be  negli¬ 
gible  in  the  released  condition. 

DESCRIPTION  OF  WING  FLUTTER  MODELS 

All  wings  tested  were  cantilever  wings  and  are  based 
on  the  section  given  in  the  following  table. 


10.  |tisf:t|l:iii»n  of  wins  I  in  vfooi  hiuh-sf tunnel.  Tin*  slop  shown  was  used  in  only  n  few  initial  test*. 
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Tin*  frequencies  of  the  vitriol  is  wings  are  given  in  the 
main  table  of  experimental  data  (table  I).  All  seetion 
constants  were  obtained  both  by  calculation  and  by  di¬ 
rect  testing.  The  basic  section  has  its  center  of  gravity 
at  42.5  percent  from  the  leading  edge.  The  stiffness 
axis  is  at  32  percent  but  was  artificially  put  at  30  per¬ 


cent  in  the  case  1  tests  by  chordwise  cuts.  (Set1  figs. 
1 1  to  13.) 

In  addition  to  obtaining  the  flutter  speed  of  the  plain 
wings,  the  effects  of  restraining  wires,  of  mass  balancing 
counterweights  in  various  locations,  and  of  large  nacelles 
both  at  the  wing  and  some  distance  away  from  it  were 
studied.  Experimental  data  are  included  in  table  I. 
In  the  aileron  tests,  the  effects  of  mass  balancing,  hinge 
location,  frequency,  and  friction  were  investigated. 

Wing  1— \V  ing  1  (see  fig.  10)  was  a  rectangular  can¬ 
tilever  wing  model  of  'j-incli  duralumin  plate  of  12- 
inch  chord  by  ',-incli  thickness  by  a  free  length  of  0 
feet  ft  inches  perforated  with  closely  drilled  h-incli 
holes  and  covered  by  a  “fnon-irich  sheet  of  duralumin 
to  give  a  smooth  surface.  The  constants  can  lie 
obtained  from  data  in  the.  experimental  table. 

Wings  2,  3,  and  4. — Wings  2,  3,  and  4  represent  a 
series  of  cantilever  wings  of  the  same  root  section  (1- 
foot  chord  by  inch  thick),  the  same  span  (0  feet  ft 
inches),  hut  having  taper  ratios,  respectively,  of  1:1, 
2:1  and  4:1.  (See  figs.  1 1,  12,  and  13.)  The  wings  are 
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made  of  (Itiraluinin  anil  are  constructed  to  give  simi¬ 
larity  in  strength  and  mass  distribution.  Note  that 
the  detail  at  or  near  the  tips  is  a  scaled-down  replica 
of  the  detail  at  the  root.  The  stiffness  axis  a  is  put  at 
40  percent  chord  or  a  -0.4  by  means  of  chord  wise 
cuts. 

The  three  types  of  wing  2  (2A,  2B,  and  2(  ’)  were  so 
designated  because  the  first  one.  2A,  finally  showed  a 


KiM  kk  II.  - Kiuliiimular  omit  ilex  ei  2\.  NoteehopJwiseeiit*  nse.|  for  pur|xw 

of  lowering  toisioiuil  frvtjwwy  -.m*l  for  pht-mg  stillness  a\is  ;»i  :«»  percent  chord 
from  lea*  ling  eilm*. 

crack  and  had  to  he  replaced  with  26,  which  is  almost 
identical.  Wing  26  finally  broke  at  the  root,  was 
repaired  by  shortening  it,  and  was  used  for  some  tests 
under  the  designation  2(\ 

Wing  5.  Wing  5  was  also  a  solid  duralumin  rec¬ 
tangular  cantilever  wing  of  1-foot  chord,  4-foot  length, 
and  I -inch  thickness  at  the  maximum  ordinate:  it  was 
used  for  aileron  testing.  (See  fig.  14)  Three  ailerons 
were  tested,  14,  24,  and  .44  inches  long  with  2,  4,  and 
4  hinges,  respectively.  Most  of  the  tests  were  per¬ 
formed  on  the  24-inch  aileron  (aileron  A  II). 

Tests  were  made  for  different  spring-restraints  on 
the  hinge,  with  a  balance  counterweight  on  the  out¬ 
board  end  (fig.  15)  and  with  a  special  arrangement 
permitting  the  changing  of  the  hinge  axis  from  the 
forward  edge  of  the  aileron  to  about  40  percent  of  the 
aileron  chord  behind  the  center  of  gravity. 


Kim  nt  12.  Tapered  cantilever  wine 
■i.  ta|MT  rat iu  hiith  in  cln.pl  ami 
thickness  is 2:1.  Dimensional  simi¬ 
larity  of  cross  section  ami  cuts. 


J  im  nt:  i:i.  Taperctl  can¬ 
tilever  ttimj  I;  tajMT 

tali".  1:1. 


Kim  kk  H.  t'anlilever  winjr  5  nse«l  for  aileron  tests. 
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Wings  6  and  7.  Wings  fi  ami  7  hit  model  wings  of 
normal  density  built  bv  covering  a  balsa  structure  with 
'(s-incli  mahogany.  Wing  6  has  the  same  external 
dimensions  as  wing  2  (fig.  10).  Wing  7  has  a  root  chord 
of  18  inches,  a  maximum  thickness  of  1.5  inches,  and  a 
taper  ratio  of  3:2  (fig.  10).  All  tapered  wings  were 
tapered  equally  in  chord  and  thickness. 

DISCUSSION  OF  EXPERIMENTAL  RESULTS 

The  scheme  already  discussed  of  introducing  flutter- 
bending  modes  completely  fits  the  experimental  results 
into  the  theoretical  picture.  Figure  17  shows  the 
theoretical  flutter  speed  for  wings  2A,  2B,  3,  and  4 
with  the  experimental  points  plotted.  Wing  2A  with 


speed ,  which  fact  again  tends  to  prove  the  contention 
that  the  flutter  bending  response  is  closely  related  to 


Fi«h  KK  Hi.  Kurt angular  wooden  wing  »'.  * I«*f1  >  anti  lit) n»r |  wing  7  right  . 

Ih|kt  ratio,  .1:2. 

the  second  ordinary  bending  mode.  Note  also  that  the 
observed  minimum  speeds  correspond  very  nearly  to 


a  flutter  speed  of  202  miles  per  hour  obviously  bends  in 
a  “first”  flutter  mode  that  approaches  the  second  bend¬ 
ing  mode  in  appearance  and  frequency  (fig.  18).  The 
flutter  frequency  calculated  on  the  basis  of  this  bending 
mode  closely  checks  the  measured  flutter  frequency 
(fig.  19).  Wing  3  checks  equally  well;  its  bending 
frequencies  are  noted  in  table  1.  Wing  4,  the  most 
tapered  one,  obviously  collapsed  (fig.  20)  in  the  second 
flutter  mode.  (See  fig.  17.)  On  this  assumption,  its 
experimental  flutter  speed  also  fits  well  in  figure  17. 

Since  the  effect  of  the  bending  mode  was  brought  so 
strongly  into  the  picture,  an  independent  study  was 
made  on  the  rectangular  wing  2B  and  on  the  tapered 
wing  3  by  attaching  one  point  of  the  torsional  axis 
rigidly  to  the  tunnel  walls  by  restraining  wires.  The 
results  are  shown  in  figures  21  and  22.  Note  that  the 
wire  attached  to  the  tip  had  no  effect  on  the  flutter 


Fioi  kk  17.--' Thcorei icul  flultiT  sjhh*<1  hasod  ou  c  instants  iiertaining  to  wings  2A.  211. 
3,  and  4.  o  =  —  0.4;  /..  =  !«;  r„*=*u.3125;  hu**=38.8  miles  per  hour.  Kxperi- 

uienta!  test  points  are  also  shown,  and  flutter  modus  and  fretpuncios  are  indicated. 


the  minimum  theoretical  speed.  (See  fig.  17.)  Of 
practical  importance  is  the  fact  that  a  stay  near  the 
root  of  a  wing  gave  a  higher  bending  frequency  and 
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¥ iii u r i  18.— Wing  2A  in  flutter,  demonstrating  first  flutter  mode.  Note  tendency  for  node  at  tip  and  maximum  amplitude  near  middle. 
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definitely  lowered  the  critical  speed.  There  existed 
points,  however,  near  the  middle  of  both  wings  for 


0  .2  .4  .6  .8  1.0  1.2  1.4  1.6 


Figure  19.—  Theoretical  flutter  frequeurtes  based  on  constants  tor  wings  2,  3,  and  4 
with  ex(>enmentully  observed  values  for  wings  2A  and  3.  *  = !  a  -  -0.4;  ia **  f «; 
r„J=0.3l25.  Case  1  (A,  a). 

which  the  stays  caused  the  flutter  speed  to  attain  a 
large  value.  The  explanation  is  that,  with  this  point 


fixed,  the  average  h  value  becomes  very  small  and 
the  h  deflection  becomes  "ineffective.”  A  relatively 
high  flutter  speed  results. 

The  matter  of  leading-edge  counterweights  has  been 
investigated,  in  particular  on  wing  2(\  Figure  2d 
shows  the  effect  of  moving  a  counterweight  along  the 
span.  The  weight  has  a  rather  surprising  iteijuiire 
effect  near  the  tip,  indicating  that,  in  this  case,  there 
must  be  an  h  node  inside  the  tip  and  again  substan¬ 
tiating  the  theory  of  the  Hotter  modes.  Farther  in 
along  the  wing  there  was  art  expected  increase  in  the 
flutter  speed.  When  all  three  weights  were  applied 
at  the  same  time,  the  flutter  speed  for  wing  2C  was 
increased  to  205  miles  per  hour,  which  is  in  good 
agreement  with  the  calculated  value. 

A  large  nacelle  at  an  inboard  position  (fig.  24)  in¬ 
creased  the  flutter  speed  from  202  to  210  miles  per  hour 
when  in  the  forward  position  and  decreased  it  to  107 
when  in  the  rearward  position. 


Figure  20.— The  cfToct  of  violent  flutter  (in  second  mode)  on  wing  4. 
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Large  bodies  disposed  at  some  distance  from  the 
wing,  such  as  floats,  were  very  detrimental  as  regards 
the  critical  speed.  (See  fig.  25.) 

Wing  1  showed  a  flutter  speed  that  is  in  agreement 
with  the  predicted  value  within  about  1  percent.  This 
agreement  is  due  to  the  considerable  internal  damping 
of  this  wing.  Wing  6,  a  rectangular  wing  of  the  same 
plan  dimensions  as  wings  1,  2,  3,  and  4,  but  of  low 
density,  showed  a  flutter  speed  about  3  percent  below 
the  theoretical  value  based  on  the  measured  parameters 
and  the  lowest  ordinary  bending  mode.  This  result 
indicates  that,  for  damped,  low-density  structures,  the 
flutter  mode  approaches  the  first  bending  mode  some¬ 
what  more  than  hitherto  indicated. 

Wing-aileron  flutter  has  been  studied  on  wing  5. 
(See  fig.  14  and  table  1A.)  The  theoretical  response  is 


shown  with  proper  constants  in  figure  20  for  the  most 
representative  aileron  All,  upon  which  most  of  tin* 
tests  were  made.  A  number  of  test  points  have  been 
directly  plotted  in  this  figure.  In  order  to  obtain 
internal  friction,  a  lead  hinge  was  used  in  some  tests. 
It  is  rather  remarkable  how  well  the  theory  is  re¬ 
affirmed  by  the  test  data.  Apparently,  if  all  param¬ 
eters  could  be  satisfactorily  determined,  no  flutter 
testing  would  be  necessary.  Tests  in  which  the  hinge 
axis  (fig.  27)  was  changed  show  the  beneficial  effect  of 
decreasing  the  uerodynamic  moment  around  the  hinge. 
The  lower  flutter  speed,  which  is  the  one  of  practical 
importance,  is  considerably  increased  as  the  hinge 
axis  is  moved  backward.  This  increase  is  not  only  a 
eenter-of-gravity  effect  but  is  also  caused  by  the  de¬ 
crease  in  the  aerodynamic  moment  around  the  hinge. 
Xote  that,  as  the  center  of  gravity  is  approached,  the 
flutter  speed  rather  suddenly  becomes  infinite. 

AIR  DAMPING  OF  FORCED  VIBRATIONS 

This  report  has  heretofore  been  concerned  with 
u  study  of  a  border  velocity  separating  stable  and 


unstable  velocity  regions.  Further  light  on  the  whole 
matter  of  flutter  is  given  by  a  study  of  the  vibration 
response  of  the  wing  system  to  impressed  forces  and 
moments,  that  is,  generalizing  the  point  of  view 


term  .\far'u'  in  c<|uutinn  (It)  a  term 

\fg>  ■’(f'tlVd/62,  and,  in  equation  ((')  a  term 

/Vi'",+*,’AV/6  are  introduced.  Here  M„  and  M,i 
are  the  magnitudes  of  the  sinusoidal  impressed  torques 
in  the  a  and  /3  degrees  of  freedom,  /'„  is  the  magnitude 
of  the  impressed  force  in  the  h  degree  of  freedom,  u  is 
the  circular  frequency  of  the  forced  vibrations,  and  the 
fs  are  certain  phase  angles. 
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from  which  both  the  nmplituiles  nml  the  plmses  iiuiv  be 
ohtainrd.  The  /f’s  ami  the  /’s  ure  listeil  in  the  appen¬ 
dix  and  A,*  and  A/  represent  the  real  and  the  imaginary 
equations  listed  under  the  calculation  scheme  for  ease  I . 

Consider  the  equation  for  A,/ and  denote  the 
determinant  in  the  numerator  by  A'»,  i.e., 

Y  H.l*  '•  ila *  lilt14" 

‘  *  -Hi.  I 


If  tlie  excitation  is  only  in  the  A  decree  of  freedom 
m  - 0,  i.  e.,  there  is  no  impressed  torque  alanit  the 
elnstir  axis.  If  a  single  exciter  were  placed,  for  example, 
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Y  i-.i  hi  ."*>  l'h**ti«cr«|>h  of  noxirl  *h*>»ine  f«e«r*|  » rr»|"«n*<*  in  !»'•  ■Ir'ctwc^ 
•  •f  frrfintti  ■lr|x>u.Jir,k‘  «>n  -iif  *|w*l  S«»lr  Outfit 

It  i«  convenient  to  define  n  certain  static  deflection 
h„,  which  is  tho  deflection  tlm*  to  force  l‘„. 


Define  a. i  ns  tin*  static  torsionnl  resjMmse  tonn  impressed 
inoliieiit  l‘J, 


Tl.cn 

•i.,  «i„’  A  »  1 

<>,,  w  A  s 

Figure  -js  is  a  photograph  of  a  three-dimensional 
moild  of  tlic  response  ratio  A  A,,  as  a  function  of  the 
cxciiinjr  frequency  ratio  u-  <*■„  ami  of  the  'peed  ratio 
r  hui.  for  a  case  of  ilellection-torsion  for  uliich  a  critical 
flutter  »|H-e<l  exists.  In  this  example,  tin*  excitint:  force 
Inis  Itecn  assumed  to  act  in  tin*  deflection  decree  of 
frcislolii  alone  i in  III.  At  /.ero  sjMi'd  there  exist  tuo 
resonant  frequencies  eorrcsfaimliiiyr  to  the  natural 
frequencies  in  the  separate  di^ree  of  freedom  The 
air  damping  due  to  speed  is  zero  ami  the  re»|>onso  is 
inlinile  at  those  freiptencie«.  (With  friction  and 
//„  the  re*|M»li*es  are  approximately  1  and  1 
respectively.!  As  the  speed  increases,  the  air  dampmy' 
increases  and  the  res|Mitise  diminishes  until,  idoin;  one 
frequency  hratieh  (the  ui„  hralich  i  the  res|Ml|se  lusonics 
nejrlttribly  small  Ahum  the  other  frequency  hratieh. 
however,  a  minimum  res|nin v  is  reached,  after  uliich 


Then 

w  here 


This  result  jlixes  the  steady-state  deflection  response 
in  t (Tins  of  the  static  deflection  A,,  due  to  an  impressed 
force  or  moment  The  results  of  some  numerical 
calculations  will  shortly  lie  presented. 

The  torsional  response  can  lie  similarly  calculated 
lad  .V.  represent  the  determinant  in  the  numerator  of 
the  expression  for 

1  mi '  *•  *>'  II „  ’•  (7,» 

|i  •  >7„  , 
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dm  air  clatii pin^  ilecrouM**  mill  llm  rospoitse  iimreiiM't*  rVitlnpln  i-ilctl.  i»  shown  in  lipuri*  .ftl.  Tin-  v nine  of 
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The  miixiinuin  response  is  t lien 


'  t_  f/fl.  ia 


Figures  35,  30,  and  37  pertain  to  titis  ease.  The 
resonant  frequency  is  strongly  affected  by  the  position 


is  obtained.  The  air  damping  increases  with  increase 
in  speed  and  the  resonant  frequency  decreases;  but  a 
speed  is  ultimately  reached  where  the  response  increases 
again  until,  at  a  vanishingly  small  resonant  frequency, 
the  response  is  very  large.  In  figure  35,  which  illus¬ 
trates  a  case  for  the  midchord  position  of  the  axis, 
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Fimri::i>v-  iVak  n*s|Minsf  A i»A.»  against  vei<nity  latio  r.Au?*  for  a  wis**  «»f  llr\»U‘*ail«Tun  flutter.  shotting  Fn.i  iti:  :«*.  --Fniml  vihialimi  t espouse  of  aileron  alone 

‘fleet  nf  ehunuing  friction  coefficients  Qh=q* = u  t-«  0a=0*=<1.2.  (Other  parameters  art*  i=0..'»:  a  =  'i.i:  against  frei|ijene>  ami  velocity ;  frielioh  coefficient 

t  j:  —  J  i  «ti ;  1.1  =  *  »«> ;  1  u\i  u: ) . 


(i  0,  tho  peak  response  occurs  at  >'>bua  —  7.10.  This 
sort  of  iiistnhility  1ms  been  culled  divergeneo.  The 
divergence  velocity  ran  he  easily  calculated  ns  follows: 
Let  1  /•  >oo  in  thi'  expression  for  a,  a,,. 

Tlii'ii 


C  ASK  2  Ifl.  M 

A  study  of  the  response  characteristics  to  forced  vi¬ 
brations  is  also  of  some  interest  in  the  case  of  flcxure- 
aileron.  The  details  are  omitted  here.  Two  figures 
are  presented.  An  illustration  showing  the  peak  re¬ 
sponse  ratio  hjli,,  in  this  case  with  and  without  friction 
is  presented  in  figure  38.  A  response  for  one  degree  of 
freedom  of  the  aileron  alone  is  shown  in  figure  3!(. 

REMARKS  ON  FLUTTER  IN  AIRPLANES 

WING  FLUTTER 

The  wing  may  flutter  as  a  whole  in  torsion-flexure. 
This  case  is  the  most  easily  treated  Kxperier.ee  with 
models  indicates  that  this  flutter  speed  may  he  calcu¬ 
lated  on  the  basis  of  the  measured  constants  with  an 
accuracy  of  a  few  percent.  The  actual  bending  fre¬ 
quency  involved  in  flutter  is  apparently  not  exactly  the 
lowest  ordinary  betiding  frequency  but  a  slightly 
higher  value. 

Probably  the  most  common  type  of  wing  flutter  is 
case  2  (flexure-aileron).  This  type,  as  well  ns  that 


involving  torsion-flexure,  is  evidently  symmetrical  with 
respect  to  the  fuselage.  The  ailerons  would  therefore 
be  in  phase  and  have  a  frequency  considerably  in  excess 
of  the  wing-bending  frequency.  This  condition  is 
favorable.  Any  slack  in  the  aileron  cables,  however, 
permits  a  motion  that,  may  cause  a  mild  type  of  flutter, 
which  should  not  be  permitted  for  too  long  a  time. 

A  nonsy m metrical  aileron  motion  would  involve  a 
second  bending  mode  (nonsymmetrical).  It  is  prob¬ 
able  that,  in  most  cases,  the  node  would  be  close  to  the 
middle  of  tho  aileron  and  therefore  poorly  coupled. 

There  remains  to  consider  a  complete  case  of  flutter 
(torsion-flexure-aileron).  Apparently  cases  do  exist  in 
which  this  type  would  appear  at  the  lowest  speed.  The 
effect  of  the  additional  degree  of  freedom  can  probably 
be  taken  care  of  by  a  safety  factor  applied  to  the  flutter 
speed  obtained  for  two  degrees  of  freedom.  The  cal¬ 
culation  of  the  case  of  three  degrees  of  freedom  is  per¬ 
fectly  straightforward  although  more  lengthy  than  the 
simple  cases. 

TAIL  FLUTTER 

lurpgard  hi  tail  flutter,  the  situation  is  more  complex. 
The  possible  combinations  are  subdivided  as  follows 
into  three  main  groups,  which  will  be  separately 
analyzed: 

(1)  Vertical  flexure  of  tail  assembly. 

(2)  Horizontal  flexure  of  tail  assembly. 

(3)  Torsion  of  tail  assembly. 

Vertical  flexure.-  It  is  possible,  in  general,  to  iden¬ 
tify  two  responses  in  vertical  flexure;  one  corresponding 
to  the  fundamental  bending  mode  of  the  fuselage  and 
the  other,  to  the  bending  mode  of  the  horizontal-fin 
arrangement.  The  frequency  of  the  fin  arrangement  is 
slightly  greater  than  the  one  obtained  with  the  rear  end 
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of  the  fuselage  fixed  in  space.  It  is  probable  that  the 
fuselage  bending  mode  need  not  be  considered.  In  any 
case,  the  flutter  speed  calculated  for  each  of  the  two  cases 
will  not  differ  very  much  because  the  density  involved 
is  about  in  the  same  ratio  as  the  squares  of  the  fre¬ 
quencies  involved.  Only  the  vertical  fundamental 
bending  frequency  of  the  horizontal  fin  will  therefore  be 
considered.  This  bending  frequency  may  couple  with 
fin  torsion  and  elevator  motion.  This  motion  is  neces¬ 
sarily  symmetrical  and  simulates  the  motion  of  the 
main  wing  system.  Since  the  elevator  has  no  particular 
restraint  to  this  motion,  it  is  evident  that  an  unbalanced 
elevator  is  highly  undesirable.  As  in  the  case  of  the 
wing  system,  the  most  frequent  cause  of  flutter  is  also 
the  flexure-aileron  combination. 

Horizontal  flexure. — Horizontal  flexure  affects  the 
vertical  fin  or  fins  and  may  be  separately  considered  as 
a  cantilever  wing  with  an  aileron.  The  rudder,  which 
takes  the  place  of  the  aileron,  has  no  particular  restraint 
to  this  motion  and  must  therefore  be  carefully  mass 
balanced. 

Torsion.- — Torsion  is  composed  of  a  relatively  low- 
frequency  type  of  flutter  involving  the  fuselage  and  a 
higher  frequency  type  involving  the  fin  tips,  which  may 
be  considered  as  fin  flexure  in  opposite  phase.  This 
type  of  flutter  is  not  common  because  the  great  stiffness 
of  the  torque  tube  prevents  its  occurrence  (0  large). 
The  rudder  is  affected  in  the  same  manner  as  for 
horizontal  flexure. 

In  summary,  it  may  be  said  that  the  expected  causes 
of  flutter  are  the  in-phase  motion  of  the  horizontal  fins 
as  flexure-elevator  and  the  motion  of  the  vertical  fin  in 
flexure  with  rudder  motions. 

With  two  rudders  disposed  at  the  ends  of  the  hori¬ 
zontal  fins,  care  must  be  taken  that  the  flexure  fre¬ 
quency  is  sufficiently  high.  The  mass  of  the  rudders 
at  the  ends  of  the  horizontal  fins  also  affects  the  param¬ 
eters  used  in  the  fin-stability  calculation;  that  is,  the 
density  of  the*  fins  ami  the  radius  of  gyration  are  in¬ 
creased.  The  eenter-of-gra vity  location  may  also  be 
changed. 

GROUND  TESTS 

DBT ICRMIX ATION  Of  CF.NTF.R-OF-ORAVITY  LOCATION 

From  the  theory,  it  may  be  observed  that  the  location 
of  the  torsional  stiffness  axis  is  of  fairly  secondary  im¬ 
portance.  The  location  of  the  eenter-of-gra  vity  axis,  on 
the  other  hand,  is  of  great  importance.  The  application 
of  a  very  low-frequency  (zero)  ton  pie  will  rotate  the 
wing  around  the  torsional  stiffness  axis  a ;  the  application 
of  a  very  high  frequency  (infinity)  torque  will  cause 
the  mass  center  line  to  remain  stationary.  As  the  tor¬ 
sional  frequency  for  wings  is  several  times  larger  than 
the  lowest  bending  frequency,  it  can  be  shown  with  nil 
desired  accuracy  that  the  axis  observed  for  the  tor¬ 
sional  frequency  is  the  center-of-gravity  line. 


Assume  the  wing  to  be  vibrating  around  an  axis  at  a 
distance  d  behind  the  stiffness  axis.  The  moment  of 
inertia  reduced  to  the  center  of  gravity  is 

I=M(ra2~xa‘)b- 

Moving  the  axis  to  d  increases  the  moment  of  inertia  to 
M(rJ—rJ)b2+  M  (xa—d)-b-  —  M(rJ  ~d:—2dxa)lr 
The  corresponding  torsional  stiffness  is  ('„■• -<p(\. 

The  frequency  is  consequently 

"  \I{rfrd2~27lx.)lr 

Then 

2  log  w  =  log  (6'«  *- d-'C,  l62— log  MlrirJ-  (P—  2<lxa) 

The  wing  will  assume  the  state  of  vibration  giving  the 
greatest  frequency.  Bv  derivation,  with  respect  to  d, 

dCh _ d  j'a  . 

(  -f-  <P(  ’ll  <P  —  2(1  j’a 

or 

d( h  21/ 

:  lO  Al 

d  —  Xa 

Then 

d-Xar  1  . 

1  —  t  Mi,  u>  r 

or,  with  d  known, 

d[l  lav,  a>c]~d  for  cc,,  '•* 

That  is,  the  center-of-gravity  axis  is  slightly  ahead  of 
the  dynamic  axis  (assuming  both  axes  to  la1  normally 
behind  the  stiffness  axis).  If  the  torsional  frequency  is 
very  large,  they  coincide  (d=ra).  (If  the  torsional  fre¬ 
quency  is  very  low.  </--(),  giving  the  stillness  axis.) 

In  other  words,  the  center-of-gravity  location  along  a 
finished  wing  can  be  determined  by  establishing  the 
dynamic  torsional  axis. 

DIMENSIONAL  CONSIDERATIONS 

Proportionally  increasing  all  dimensions  of  a  wing 
while  retaining  all  details  lowers  the  frequencies  in  in¬ 
verse  proportion  to  tin*  size.  The  reference  speed  ub 
therefore  always  remains  the  same,  as  do  all  other  pa¬ 
rameters  including  the  wing  density.  The  actual  flutter 
speed  therefore  depends  on  the  shape  but  not  on  the  size. 
It  is  important  to  keep  in  mind,  however,  that  the 
reference  is  to  wings  or  tails  similar  in  all  respects.  In 
reality,  a  lighter  eon  t  met  ion  is  necessarily  employed  in 
larger  wing  sizes,  resulting  in  a  weaker  structure  and  a 
general  lowering  of  the  critical  flutter  speed. 

The  foregoing  considerations  are  significant  in  the 
testing  of  models.  Thus  a  true  model  constructed  of 
the  same  material  as  the  full-scale  airplane  will  have 
the  snmc  flutter  speed.  For  testing  purposes,  it  is 
very  desirable  to  have  a  fairly  low  flutter  speed.  This 
end  may  be  achieved  by  employing  models  of  special 
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unit  (‘rials  related  to  celluloid,  which  have  a  value  of 
V/vp  nearly  five  times  smaller  than  that  of  materials 
normally  used  in  airplane  construction. 

The  most  desirable  condition  would  he  to  use  a 
material  with  the  same  density  p  as  the  airplane  and 
with  the  moduli  K  and  G,  say,  1/n  times  the  original 
values.  Fortunately,  the  density  of  the  model  wing 
can  be  very  simply  corrected  by  using  a  suitable  thick¬ 
ness  of  the  materials.  Thus,  if  the  density  of  the 
material  used  is  three  times  lower  than  that  of  Un¬ 
original.  the  thickness  of  the  skin  and  all  the  internal 
members  is  increased  by  a  factor  of  3. 

It  should  be  further  noted  that  the  model  can  be 
critically  cheeked  as  to  accuracy  of  reproduction  In¬ 
direct  measurements  of  its  mechanical  properties.  In 
other  words,  all  the  parameters,  including  the  reference 
quantity  ub.  are  directly  measured  on  the  model  itself. 
The  value  of  ub  is  usually  close  to  the  predicted  value. 
The  important  point,  however,  is  that  it  is  not  neces¬ 
sary  to  depend  on  a  predicted  theoretical  value. 

Thus  the  feasibility  of  conducting  direct  flutter  tests 
on  models  of  actual  airplanes  or  of  its  component  parts 
is  indicated.  Some  work  of  this  nature  is  now  being 
undertaken.  The  procedure  may  be  of  value  in  cases 
that  are  difficult  to  treat  theoretically  and  should  be 
of  value  in  accumulating  useful  experience  on  special 
designs. 

GENERAL  CONCLUSIONS 

1.  The  two-dimensional  theory  has  been  verified 
within  the  limits  of  error  in  the  determination  of  the 
primary  parameters. 

2.  The  most,  essential  three-dimensional  effect  is  (lie 
occurrence  of  distinct  flutter  bonding  modes,  which 
differ  from  the  ordinary  vibration  modes  in  that  they 
tend  to  assume  a  form  which  approaches  the  next 
higher  vibration  mode  and  exhibit  a  correspondingly 
higher  frequency.  Tin-  flutter  speed  is  consequently 


lower  than  that  calculated  on  tin*  basis  of  the  lowest 
vibration  frequency  and  the  flutter  frequency  itself  is 
higher.  For  ordinary  damped  structures,  this  effect 
lowers  the  flutter  speed  calculated  on  the  basis  of  the 
lowest  bending  mode  by  only  a  few  percent. 

3.  A  cantilever  wing  flutters  at  a  speed  calculated 
by  using  the  constants  for  the  most  representative 
section,  which  is  located  at  approximately  three- 
quarters  of  the  semispan. 

4.  Aspect  ratio  and  structural  damping  effects  tend 
to  increase  the  flutter  speed  bv  a  few  percent  above 
that  calculated  for  infinite  aspect  ratio  and  zero  internal 
damping. 

">.  The  effect  of  mass  balancing  to  bring  the  center 
of  gravity  forward  is  essentially  as  predicted  by  theory. 
The  effect  of  nacelles  is  of  lesser  importance,  but  large 
weights  loeated  at  some  distance  away  from  the  wing 
and  attached  to  it  show  a  very  detrimental  effect  on 
the  flutter  speed. 

(I.  Wing-aileron  experimental  studies  show  that  the 
characteristic  flutter  range  predicted  by  the  theory 
exists  and  is  in  substantial  agreement  with  the  pre¬ 
dicted  values.  A  decrease  in  the  unbalance  and  an 
increase  in  the  frequency  ratio  are  both  beneficial. 
There  exists,  for  each  value  of  the  unbalance,  a  certain 
critical  frequency  ratio  and,  inversely,  for  each  fre¬ 
quency  ratio,  a  certain  critical  value  of  the  unbalance. 

7.  The  considerable  difficulty  involved  in  the  deter¬ 
mination  of  the  primary  structural  parameters  includ¬ 
ing  the  damping  is  recognized  and  will  constitute  one 
of  the  chief  problems  of  future  flutter  research. 


Langley  Memorial  Aeronautical  Laroratory, 
National  Advisory  Committee  for  Aeronautic^- 
Langley  Field,  Va.,  St/itanbtr  22.  I  OSS. 


APPENDIX 


LIST  OF  NOTATION 

a,  angle  of  attack  (fig.  40). 

/S,  aileron  angle  (fig.  40). 

h,  vertical  distance  (fig.  40). 

da  ..  (Pa  , 
a==  Wa^dP,ctc- 

a0,  /So,  h0,  <fiU  if amplitudes  and  phase  angles. 

b,  half  chord,  used  as  reference  unit  length. 

a,  coordinate  of  axis  of  rotation  (torsional  axis) 

(fig.  40).  location  of  stiffness  axis  in  percentage 

total  chord  measured  from  the  leading  edge,  is 

... _ 2  (stiffness  axis)  , 

100  2  or  a= - - 1. 

c,  coordinate  of  aileron  hinge  axis  (fig.  40).  Location 
of  aileron  hinge  axis  in  percentage  total  chord 

1  +c 

measured  from  leading  edge  is  100——  or 

_2  (aileron  hinge) 
r  100  • 
p,  mass  of  air  per  unit  of  volume. 

Af,  mass  of  wing  per  unit  span  length. 

7 r 

«=  ~\j~>  the  ratio  of  the  mass  of  a  cylinder  of  air  of  a 

diameter  equal  to  the  chord  of  the  wing  to 
the  mass  of  the  wing,  both  taken  for  equal 
length  along  the  span;  this  ratio  may  be 
expressed  as  k=0.24(62/U’)  (p/p0)  where  H  is 
weight  in  pounds  per  foot  span,  b  is  in  feet, 
and  p/po  is  ratio  of  air  density  to  standard 
air.  [The  quantity  a—  H7461  (weight  per 
square  foot  per  chord  in  feet)  lias  been  used 
by  British  writers.  Thus,  k  (<l.0(i/<r)  (p'p„).l 

/„=  j ^ •  location  of  center  of  gravity  of  wing-aileron 
system  meaxurcd  from  a  (fig.  40);  S„,  static 
moment  of  wing-aileron  per  unit  span  length 
referred  to  a.  Location  of  center  of  gravity 
in  percentage  total  chord  measured  from  the 

leading  edge  is  1001  —  or  u-j xa~ 

2  (center  of  gravity)  _ 

'  100 

fj—  reduced  location  of  center  of  gravity  of  aileron 
referred  to  c  (fig,  40)'.  S$,  static  moment  of 

aileron  per  unit  span  length  referred  to  c.:l 

I  / 

r„  =  ^  radius  of  gyration  of  wing  aileron  referred 
to  n  (fig.  40)  ’  /„,  moment  of  inertia  of 

wing  aileron/about  the  elastic  axis  per  unit 
span  length. 

re~  J  ifp’  r'd  need  radius  of  gyration  of  aileron  referred 
to  c  (fig.  40).  I9,  moment  of  inertia  of 

aileron  about  c  per  unit  span  length. 

1  Note  that  M  refer*  to  the  total  wing  mass  and  not  to  the  mass  of  the  aileron  alone. 


Ca,  torsional  stiffness  of  wing  around  a  per  unit  span 
length. 

Ca,  torsional  stiffness  of  aileron  around  c  per  unit 
span  length. 

Ca,  stiffness  of  wing  in  deflection  per  unit  span  length. 

w„  natural  angular  frequency  of  torsional  vibra- 

v  “  tion  around  «  in  vacuum  (&>„  -  ‘2rfa,  where 
fa  is  in  cycles  per  second). 

I'C 

u-.i=  w  natural  angular  fmpieney  of  torsional  vibru- 
’  3  tions  of  aileron  about  c. 


ull,~  V  T7’  na*llra^  frequency  of  wing  in  detler- 

’  ‘  tioB. 

*****  Trs%7 

-/  -1/2  0 - c - -  / 


Axis  of  rotation 


1^  ^  I  e  g.  of  entire  wing 

eg  of  aileron 

Ficukk  40.  -  Half  chord  b  is  used  jls  tin-  unit  length.  Tin*  positiu-  directions  « f  <*. 
0,  and  h  are  ind United  hy  arrows.  Note  that  «  is  measured  from  uiidehord  and 
that  .r„  is  measured  from  the  clastic  axis  positive  to  the  right.  Also  note  t!mf 
is  a  “reduced"  parameter  and  not  the  actual  distance  from  the  hinge  to  the  c.  tr.  of 
the  aileron. 

(,  time. 

r,  speed  of  forward  motion. 
rr,  flutter  or  critical  speed, 
u.xircuhtr  frequency  of  wing  vibrations. 

k=  ,  ->  reduced  frequency  =  number  of  waves  in  the 

wake  in  a  distance  equal  to  the  semichord X2*\ 

\/k,  reduced  wave  length  =  length  of  one  wave  of  the 
wake  in  terms  of  a  distance  equal  to  the  semi- 
chord  X2ir. 

/'’and  G,  functions  of  k  in  table  2. 

R..=- 

R.,=  -•'«  +p.l..+(|+«  )(j,'  ~  t u) 

Rah=-Aai- l"^2^’a^r 
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1 th  ■- 

Tlic  quantities  A„,,  A*,,  etc.  and  ft„A',  R».Y,  etc.  arc 
defined  under  the  calculation  scheme  (pp.  5-7).  The 
T’s  are  listed  in  table  2.  The  definitions  of  the  Ts 
are  given  in  reference  1,  page  5,  should  other  values 
than  those  listed  in  the  table  be  required. 
vjbuu,  flutter-speed  coefficient  (cases  1  and  3). 


v/buk,  flutter-speed  coefficient  (ease  2). 

<jijua,  frequency  ratio  (case  1). 

aifi/uh,  frequency  ratio  (case  2). 

up/uu,  frequency  ratio  (case  3). 


1/2 

(1/2) +«’ 


divergence  velocity. 


9a.  9a.  9>„  structural  damping  coefficients;  *9  corre¬ 
sponds  approximately  to  the  usual  logarithmic 
decrement. 

Ma,  M3,  magnitude  of  sinusoidal  impressed  torques  in 
the  a  and  &  degrees  of  freedom. 

/V  magnitude  of  impressed  force  in  the  h  degree  of 
freedom. 

ip,  phase  angle. 

a0'asl,  fi0 h0  7i,„  peak  response  for  the  various 
degrees  of  freedom. 
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Values  of 
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-0.3 

-0.2 

-0.1 

0 

0 

1.00000 

0.  90250 

0.81000 

0. 04000  i 

0.  490(H) 

0. 30000  ! 

0.  2.5000 

.  100 

!  i.ooooo 

.  90256 

.81014 

.64031  | 

.490  U) 

. 36070  1 

.  25094 

. 10007 

1.00000 

.90268 

.  81038 

.  64083  i 

.  49135 

.36192 

.  25258 

.250 

1.00000 

.90291 

.  81086 

.64187  1 

.  49303 

. 36432 

.  25578 

.  33333 

1.00000 

.90323 

.  81154 

.  64331 

.  4953.5 

.  36764 

. 20019 

.  500 

1.00000 

.90415 

.81345 

.64741  1 

.  50189 

.  37089 

.  27240 

.  06667 

1.00000 

.90524 

.81570 

.  6522 1 

. 50902 

.  38784 

.  28690 

.  83333 

1. 00000 

.00690 

.  81922 

.65902  1 

.  52120 

.  40404  . 

.  30810 

1.000 

1.  <10000 

.  90879 

.  82297 

.66751  i 

.  53363 

.42133 

.  33060 

!  1.250 

1.00000 

.91192 

.  82938 

.  08093  : 

.  55404 

.4.5051  : 

.  30*55 

1  1.51510 

1.00000 

.91568 

.  83700 

.09691  1 

I  . 57950 

.48501 

.  40325 

1.001567 

1.00000 

.  91797 

.84175 

.  TIKKU 

.  5!Mf.!l 

j  . 50590  , 

.  44025 

1.  7*572 

1.00000 

.  91986 

. 84557 

.71457 

1  .00701 

.  52280 

.46215 

2.  000 

I  1.00000 

.92334 

.  85266 

. 72923 

.  02972 

i  .55413 

.  50245 

2. 27273 

1.  00000 

.92791 

. 86202 

.  74856 

.  05902 

,  .  59520 

.  5.5530 

2.500 

1.00000 

.  93187 

.  87000 

.  76500 

.  08.500 

!  .  63000 

.60000 

2.  94118 

1.00000 

.93964 

.  88576 

.  79739 

.73190 

.  69829 

.  68755 

3.  33333 

1.000(4) 

.94657 

.  89980 

.  82020 

. 77920 

.  75960 

.  76500 

4.  16607 

1.00000 

.96120 

.  92949 

.  88075  ! 

.  87206 

.  8853.3 

.  92655 

5.000 

1  1.00000 

.  97518 

.  95763 

.  94437 

.  90021 

1.00515 

1.07920 

6.250 

1.00000 

.  99478 

.  99720 

1.02490 

1.08312 

1.  17180 

1.  291 10 

8.  33333 

1.00000 

1.0235 

1.05,507 

1.  24238 

1.20195 

1.41378 

1.  59785 

10.  000 

1.00000 

1.0433 

1.09484 

.  1.32296  • 

1.38430 

1.  57904 

1.80700 

12.  500 

1.00000 

|  1.0688 

1.  146)3 

1.  42667 

1  54103 

,  1.79101 

2.  07480 

16. 66667 

1.00000 

i  1.  1029 

1.  21474 

1.56510 

i.  rfli2« 

2. 07304  1 

2. 43050 

'  20.000 

1.000(10 

1.  1246 

1.  25838 

1.05312 

1.  88422 

j  2. 25168  < 

2.  655.50 

TABLE  -I 
Values  of 

^  775*=  «dl  4  «M" 


0 

0.  1 

0.2 

0.3 

0.4 

0.5 

U.  6 

0.7 

0.8 

0  9 

1.0 

. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  100 

.(H)  123 

.00100 

.00078 

.00060 

.  00045 

.  IHKW2 

.  (KHI21 

.  (HH11 1 

.  (XM105 

.00001 

0 

. 16667 

. 00342 

.  00267 

.00219 

.00168 

.00124 

.  (HH)87 

.  IK  1057 

.IHH132 

.  (HKIl  4 

.IMMMI4  0 

.  250 

.(H)768 

. 00623 

.  (HM93 

.  00378 

.(H1279 

.IKU96 

.iH.ia; 

.(HHI71 

.(HH)32 

,  IN  MM  IS  0 

.  33333 

.01367 

.01109 

.00877 

.  00673 

.00498 

.  (M1349 

.  00225 

.  001 19 

.  000.50 

.  (MM 05  0 

.  500 

.  03085 

.02505 

.01978 

.01519 

.01122 

.  (X)787 

,(H)507 

. 00290 

.(MU ‘26 

.  (HHI33  0 

.  66667 

.  05494 

.04461 

.  03525 

.  02706 

.01999 

.01402 

.  0090.5 

.0051 .5 

.  0022.5 

.  (MM >59 

0 

.  83:133 

.08647 

.07018 

. U554U 

.04250 

.03137 

.  02236 

.01419 

. 00*06 

. 00352 

. 00092 

0 

1. 000 

.  12517 

.  101.58 

.08015 

.06146 

.  04535 

.0:1177 

.« ravi* 

.01164 

.09510 

.  (HM33 

l> 

1 . 250 

.  19729 

.  mm 

. 12621 

.  09670 

.07130 

.  0491*0 

.03214 

.01*25 

.  (HI798 

.(K)20* 

0 

1.51516 

.  29270 

. 23723 

. 18699 

. 14317 

. 10547 

. 07373 

.04744 

. 02690 

.01175 

.00507 

0 

1.66667 

.35614 

.288.54 

.  22733 

. 17398 

.  12810 

.  08952 

.  05755 

.  0.3262 

.01425 

.  (M)37 1 

0 

1.78572 

.41062 

. 33258 

.26195 

.  2<H)40 

.  147.50 

. 10302 

. 06620 

.037 50 

.  01637 

.  (KI426 

0 

2.  OOO 

.51910 

. 42022 

.  33078 

. 25290 

.  18595 

. 12980 

.  08335 

.04717 

.020.5* 

.(M>535 

0 

2. 27273 

. 67692 

.  54753 

. 43071 

. 32902 

.24177 

. 16857 

.  10813 

.  (Hi 1 1 3 

.  02*164 

. 00691 

0 

2.  500 

.  825.50 

. 66738 

.  52463 

.  4(H)50 

.  29407 

.2(H*8 

.  13131 

1174 1 7 

.  03230 

. 00*38 

0 

2.94118 

1.  15907 

. 93605 

. 73500 

.  56040 

.41094 

.  28.588 

. 1*296 

.  10318 

.  044K6 

.01162 

0 

3.  :im3 

1.  50670 

1.21566 

. 95364 

. 72637 

.  53204 

.  *>968 

.  236*) 

.  13308 

.05779 

. 01495 

0 

4. 16667 

2.  40775 

1.93923 

1.51852 

1.  15435 

.  84375 

. 58495 

. 37303 

.  20958 

.09080 

.02342 

0 

5.  OOO 

3.  53385 

2.  84190 

2.22193 

1 . 68628 

1. 2*136 

.  85134 

.54187 

.  *t:i8o 

. 13137 

.03:4*0 

0 

6.250 

5.  65191 

4.  (HW26 

3. 54031 

2. 68138 

1  95216 

1.34764 

.  8.5572 

.  478.5:1 

.2IH144 

. 05296 

0 

8.  33333 

10.  34067 

H.  28124 

6. 44720 

4.  87081 

3.  53670 

2.  43454 

1  .54135 

.85928 

. 36900 

.09447 

0 

10.  000 

15. 14353 

12.  11112 

9.41589 

7. 10321 

5.  14958 

3.  53887 

2. 23667 

1.24464 

.53443 

.13633  1  0 

1 

TABLE  4— Continued 
Values  of  /«" 

I M  1,1/  ha" 


\  c 
l/k 

0 

0. 1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.8 

0.9 

1.0 

0 

ft.  32298 

0.  23827 

0.  169*) 

0.11472 

0. 07306 

0. 04275 

0. 1)2213 

0.  00944 

0. 00283 

0.  00035 

.  100 

.  32289 

.  23820 

. 16925 

. 11468 

.07304 

.  04273 

.02212 

.  00944 

.IMJ2S3 

.  (HH)35 

0 

. 16667 

.  32273 

.  23808 

.  16916 

.11462 

.  07*H) 

.  04269 

.02211 

.  (H1943 

.  (H12H.4 

.  (XM)35 

0 

.250 

.  32242 

.  237K4 

.  16899 

. 11450 

.  07292 

.  04260 

.02209 

.  IHW42 

.  IM12K2 

.  (HHK45 

.  xm 

.32198 

.  23751 

.  16874 

.11433 

. 07280 

.  04243 

.  02205 

,  009W 

. 002*2 

.  00035 

.  5(H) 

. 32075 

.  23658 

.  16806 

.11:486 

. 07250 

.04241 

.02195 

.  IKWI7 

.  IH)2S0 

.  (HHK45 

. 66667 

.31931 

.23548 

.  16727 

. 11331 

.07213 

.  04220 

.02184 

.  00931 

.  IH)279 

.  (HH):14 

.  *4333 

.31714 

.23385 

.  16609 

. 1 1249 

.07161 

.04183 

.  02168 

.  <K>925 

.  00277 

.  (MMI34 

1  (HH) 

.31483 

.2:4212 

. 16484 

.11163 

.07106 

.04156 

.02150 

.  (H1917 

.  (H127I 

.(MM)34 

1.250 

.31090 

,22919 

. 16273 

.11019 

.07013 

.04101 

. 02121 

.  (H)9(>4 

.  (HI271 

.  (KH134 

0 

1.51516 

. 30625 

.22572 

. 16024 

.  10849 

.  06903 

.040:46 

.02087 

.00890 

.  00267 

.  (HHW3 

0 

1.66667 

. 30343 

.22362 

.  1.5873 

.  10746 

.  06837 

.03997 

.  02(N>7 

,  (HIH8I 

.  (Hrjm 

.  (KKW3 

0 

1.78572 

. 30112 

.22190 

. 15752 

.  10662 

.  06783 

.03966 

.  020.50 

.  (MI874 

.  (H)262 

.  mm 

0 

2.  (HH) 

.29688 

.21876 

. 15526 

.  10509 

.06685 

.  (13908 

.  02021 

.00861 

.  00257 

.00033 

2. 27273 

.291*) 

.21462 

.  15231 

.  10307 

. 06556 

.03832 

.01981 

.00845 

.00252 

.  (HHW2 

0 

2.5(H) 

.28657 

.21112 

.14980 

.10137 

.06448 

.014768 

.01948 

,  00831 

.  00249 

.(H1031 

0 

2.  94118 

.27726 

. 20423 

.14491 

.09803 

.06235 

.03643 

.01883 

.00803 

. 00240 

.00030 

0 

3. 

.26899 

. 19813 

. 14056 

.09,509 

.(WHO 

.03533 

.  01826 

. 00778 

.002142 

.  mrM 

0 

4. 16667 

.25166 

.  18536 

.13152 

.08894 

. 06655 

.03305 

.01708 

.00728 

.00217 

.  00027 

0 

5.  (HH) 

.23522 

.17325 

.12292 

.08315 

.05286 

.03089 

.01590 

.00681 

.  00203 

.  00025 

6.  250 

.21229 

.15642 

.11100 

.07510 

.04775 

.02791 

.01442 

.00615 

.  0018.3 

. 00023 

0 

8.33333 

.17895 

. 13197 

.09372 

.06345 

.04038 

.02362 

.01220 

.00521 

.00155 

.(H»19 

0 

10.000 

.15613 

.11525 

.08194 

.05551 

.03535 

.02069 

.01070 

.  00457 

.00135 

.00017 

0 

I 


\ 
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c.  g.  location!  fe  J  v 
■rceut  ~  ~ 

(1  froii 
.  K.) 


r. 


I  <|*C  rceut 

chord  from 


Cx 

k  i 


2  5 

w  * 

~  c 

3 

I  1 


t  1 1  S 


I.Otf  ,  «.»  19.7 


41.2  41.2  30.0  0.  3125 


41.2  12.3  .  30  n  0.3125  30.50 

i 


11*.  7 
17. 0 


11  2H 


3S.  US 


Test  conditions  and  remarks 


'I'csts  on  wing  i  were  preliminary. 
A  Stop  was  Used. 

1  Later  tests  employ ed  a  screen  instead 
Damping  test 
-  Hun  20-3  repeated 
;  Weight,  u.oan  ih.  on  L.  K.  is  in.  from 
tip. 

Stop  used;  all  other  tests  screen  used 
(tip  flutter;. 

Model  supiHirt  fractured 
Model  repaired 

Best  raining  wires  17  in.  from  tip. 
3.0  in.  from  L.  E. 

Nacelle  C.10  lb.  placed  2  ft.  Irom  top. 

2.70  in.  ahead  of  c.  g. 

Nacelle  placed  2  ft.  from  top,  1.34  in. 

back  of  c.  g. 

Hun  20-10  repealed 
Counterweight  3.07  lb.  at  L.  K.  13  in. 
from  tip. 

Hun  20-20  rejteated 

Plain;  no  weight  attached 

Disk  12.77  Ih.  12-in.  diam.  on  bottom: 

c.  g.  of  disk  on  c.  g.  of  airfoil. 

Disk  12.4.1  lb.  17-in.  diant.  in  same  po¬ 
sition  as  20-24. 

Disk  as  in  2n-2.1  placed  1  in.  ahead  "f 


-  ■ 

. _  . 

2o  2 

2.2*3X  10- 

97.  1 

20-3 

2 

99.  2 

21 M 

20-5 

2.  274 

9s.  1 

20  S 

2.  201 

111.  4 

20  10 

2. 250 

15  0 

iss.o 

20  11 

203.0 

20  1.1 

2. 154 

10.2 

202.0 

20-10 

2.210 

151.  0 

20-17 

2.  110 

9.  45 

2  !0.  o 

20  is 

2.  170 

197.  0 

:  20  19 

2.  225 

151.0 

20  20 

2.  173 

197.  o 

:  20  21 

2. 174 

2<»2  o 

1  20-23 

2.  159 

10.  2 

201.  \t 

1  20  21 

2.204 

s  31 

111  2 

20  21 

2.  22s 

7  31 

121  1 

1  '  •  K. 

i  Disk  as  in  20-25  placed  1  in.  hark  of 
'  e.  g. 

I  Hun  20-23  re  |  *eated 
Airfoil  tip  restrained  by  wires  , 

Hun  20  25 rc|*catcd 
:  Him  20-23  re|H‘itled 
Same  as  20-31.  movies  taken  . 

;  Plain;  restraining  w  ires  s  in.  from  tip 
i  Plain;  restraining  w  ires oo in. from  tip. 

Hestraining  wires  40  in.  from  tip; 

!  model  broke  near  tip. 


20  JO  2.  232 

20  27  2  211 

Jn  Js  2.  I  to 
20-20  2.  I.'.! 
2o  ;io  2.232 

20-31  2  140 

2032  2144 
2o  33  2  122 
20  31  2.  171 
2u  3.1  2  013 


lo  s 
It.  1 
s  21 


IP.  2 
131  9 

2*1.1  «i 

201  I 
120  7 
P«  .{ 
P.*s  7 
20s  o 
I  V*  i 
2M  7 


Khmer  velueitx  >m. 


in  n  . 

2.30 

12.  S 

31.2 

19.2 

42.  5 

45.  2 

1C  7 

4. 112 

25.  0 

67.2 

,  35. 8 

1  i  ! 

42  5 

43  9 

i us  of  pyrnt 
notation) 


OF  FLUTTER 


4.") 


I\S  OF  FM'TTKH  I X  V  KSTIO  AT  ION  ( 'out  inu.-.l 


l 


Tost  comliiions  and  remarks 


Model  2H  repaired 

20-73 

2.133X19- 

220  { 

\\  ilh  3.0.5-lb .  ceiunterwciphi  1  in.  from 
tip  at  Jo  K. 

20-74 

2.  169 

196.  * 

(’ounUTweipht  13  in.  freun  tip  at  L.  K 

218  75 

2.  152 

!  am. 

.  -  . . 

('ounteru  eight  25 in.  from  lip  at  I..  K 

218-76 

2.  lot 

214,  e, 

Counterweight  37  in.  from  tip  at  I..  K 

20-77 

2.  1 1 9 

215.  7 

Counterweight  49  in.  fnmi  tip  at  1..  K 

20-7* 

2.  1 22 

2*2.  2 

Coimterweipht  61  in.  from  tip  a?  I..  K 

20-70 

2.  136 

219.  * 

With  weights  simultaneously  at  13  in- 
37  in.,  and  til  in.  from  tip. 

20  M 

2.  056 

291  3 

—  —  — . 

— 

—  . 

- 

20.  48 

Plain  tapered  w  inp 

20  22 

2.  !0n 

16.4  216.0 

Kun  20-22  repeated  . 

20-iVs 

2.  150 

211  9 

Nacelle  ti.io  lb.;  c.  if.  3.1(1  in.  back  of 
v.  p.  of  airfoil. 

20-69 

2. 212 

109.  2 

Nacelle  plaecel  1.10  in.  l>ack  of  c.  p. 
Kip.  25. 

20-70 

2  245 

91.0 

Nacelle  placed  2.80  in.  ahead  of  c.  p 

20  71 

2  240 

*  .VI  90  1 

Kun  25  68  repeated 

20  72 

2  149 

14  0  213.3 

16. 8h 

Violent  flutter  suddenly.  Whip  lent 
out  of  shape.  Kip.  2t"). 

20  67 

2.  171 

203.  9 

20.48 

Restraining  wires  59  in.  from  tip  at 
stillness  axis. 

20  v* 

2  1*7 

211  1 

Same  except  wires  2  in.  ahead  of  stitf. 
ness  axis. 

2o  9o 

2.  200 

204.4 

Same  except  wires  4  in.  hick  of  still¬ 
ness  axis. 

20  91 

Restraining  wires  40  in.  freon  tip  at 
stillness  axis. 

20  93 

2  214 

155.8 

Kcstraininp  wires  17  in.  from  tip  at 
stillness  axis. 

20  91 

2  229 

152.9 

Kcstraininp  wires  I  in.  from  tip  at 
stillness  axis 

20  95 

2  182 

211.4 

Kcstraininp  wires  27  in.  from  tip  at 
stillness  axis. 

21)  96 

2  202 

200.8 

3.  45 

Kcetanputur  wooden  model;  ■d  in. 
span. 

20  S2 

2  26s 

73  8 

6.01 

Tap-reel  wooden  model;  completely 
ijcslroycd  by  violent  flutte  r. 

20  92 

2.  199 

205. 6 

itn  \'}).)o|,ia  J  »uni.4 


JNCLASSIFIED 
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Wing 

model 


Aileron  (see  fig.  14) 


|  First  'Aileron 
Wing  Ailerou  **nt-  ire-  (t>f. 


weight  weight 
(lh.)  (Ih.) 


ing  uuency  . 
(ey-  (cy* 
eke*/  ties/  llwnJ 
see.)  sec.) 


All . 

-  44.  M  0. 9So 

JO  6*  0 

AV .... 

.  m» 

0 

AH  _ 

. yso 

5.  75 

AV... 

.61* 

7.91 

All 

.yso 

12.  5 

All 

K  50 

All 

II. U 

All 

II.  II 

All. 

12.  5 

All 

13.  1 

All 

9.07 

All 

9  17 

All 

All 

II 

All 

13  3 

AH 

III. » 

All 

0 

All 

0 

All 

0 

All 

II 

All 

3.  on 

All 

6.92 

A I 

567 

1) 

AMI 

3INI 

ll 

AH! 

1.43 

II 

A 1 II 

0 

Al  11 

<1 

All! 

(1 

A  111 

4  67 

A 1 1 1 

6.0H 

.win. 

2. 01 

(1 

.Will 

II 

A  VIII 

II 

Will 

II 

A  VUI.. 

7  00 

Te<l  coin  til  Urns  and  remarks 

Aileron  with  3  free  hinges  . 

Aileron  with  springs  . 

Theoretical  results  for  uileron  AH  shown  in  fig.  26 

Hinges  free . 


3  free  hinges  with  uileron  counterweight.  o.55 
Ih.,  Fig.  IS. 

Counterweight,  0.452  Ih.  . 

Counterweight.  0.346  Ih . . 

Counterweight.  i>.3Hh  Ih 
do  ... 
do . . 

2  free  hinges. . 

3  free  hinges 

4  free  hinges 

4  free  hinges  with  le-nl  for  damping 
do  . 

...  do 

.  do 
do 

Hinge-pin  |>osition  varied  as  shown  in  fig.  27 
do 
do 

.  do  . 

...  do 


I^ower 

(Til  ic'd 

i.  Air  dens*  velocity  Air  denv 

'  lly.e  (III.  p.  Itv.#, 

It.  I 


iTItlCiiJ 

\  el«*c»t  > 
mi.  p. 
Ii.i 


a>  ;jo 

2  27:4X10-4 

5i  9 

•JO  37 

2.  27k 

5ii  7 

3h  3* 

2.  273 

Ui.  4 

JO  39 

2  271 

42.4 

JO  40 

2  3#9 

76.  h 

at  41 

2  306 

44. K 

JO  42 

2.  3Ul 

55.  4 

3f  43 

2  240  X  HI 

147  o 

a»  44 

2  229 

1 51  6 

•JIMS 

2.259 

9K  4 

JHd 

2  2M 

57.  4 

2.2:4#i 

1 5y.  .4 

a>  47 

2  26* 

5 1  II 

2  162 

19.4  o 

JO- 4* 

2  245 

52  6 

2.  165 

199.9 

ai-49 

2.  242 

5k  3 

2  US 

225  1 

acsu 

2  242 

011  3 

2. 3*n 

12-5  7 

34-51 

2  240 

65  0 

2  199 

1 49  '.# 

JO  52 

No  flnuer 

au-ftt 

\i»  dulter 

a>  54 

2.  249 

75  2 

2.  243 

39*  1 

a»  -.v. 

2  243 

5 

2-  173 

197  k 

ai-.vj 

2  243 

«*  9 

2  19.1 

167  .4 

at  57 

2  262 

6J  k 

2  21S 

147.  1 

34  5k 

2  312 

6*«  3 

2  164 

21*  7 

34  Sy 

2.  259 

57  k 

2  221 

In.’.  7 

34-6ci 

2  J5n 

#41.  2 

2  161 

3*5.  i 

3>  6| 

2  25 1 

2 

2-  199 

145  i 

Jv  6.1 

2.  3d 

7k  it 

2.  171 

!•#■#.  5 

34-64 

2.  2-V* 

97.  k 

2  Ik* 

22.4  2 

34-65 

2  259 

67  3 

2  I7k 

1*<6  K 

34  66 

2.  3Vn 

64  k 

2  21  o 

l#*|  6 

3»  S3 

2.  2*0 

64  4 

J  17k 

210  5 

34-K4 

2  JkI 

91  (. 

2  IH 

3«6  o 

34  K6 

Ni.  flutter 

34  *7 

2  2*» 

K5  5 

2  i#9K 

290  K 

34  ** 

2.  259 

4t  1 

2  l.to 

2**4  1 

j 


Positive  direction*  of  mc*  and  annl<'»  (force*  am)  moment*)  am  aliown  by  arrow* 


[  Axl* 

Pome 
ffiarallal 
to  axis) 
aymtjol 

Moment  n t mu t  anh 

■ 

Angle 

Volocltli* 

1  ><:*lgnatloil 

Hy  iii- 
boi 

Designation 

.Sym¬ 

bol 

Positive 

direction 

Doafgria- 

lion 

Hym- 

bol 

Linear 
(compo¬ 
nent  along 
axia) 

! 

Angular 

Longitudinal . 

,V 

.V 

Rolling . 

h 

Y  , 7. 

Itotl . 

<• 

u 

V 

1  literal. . 

>' 

Y 

Pitching.... 

M 

z, . *X 

Pitch.... 

0 

f 

1 

Normal .  . 

/ 

Y 

yawing.... 

X 

X  ,Y 

Yaw . 

* 

IP 

T 

Absolute  rnefficientH  of  moment 
,,  A  ,,  M 

‘  >fiS  C*Mqc8 

(rolling)  (pitching) 


(yawing) 


Anglo  of  set  of  control  surface  (relative  to  neutral 
position),  A.  (Indicate  surface  by  proper  sulwmripl.) 


4.  PROPELLER  HYMBOL8 


IK 

1  liametor 

IK 

(ii  omolric  pitch 

vUK 

I 'it'll  ratio 

C. 

V\ 

Inflow  velocity 

V., 

Slipstream  velocity 

n, 

T, 

Q, 

T 

Thrust,  absolute  coefficient. 

Torque,  alwolute  coefficient 

n, 

*, 

Power,  alrsoluto  coefTieient 

•Speed-power  coefficients  yj 
Efficiency 

Hcvolutiona  per  second,  r.p.a. 
Effective  helix  angle  «>  lair 


I  bp  7f‘t,04  kg*m/a»»5C0  fHh./acc. 
I  metric  horsepower'*-*  I  .0132  hp. 

1  iri  p.b,  0,4470  in.p.i«. 

1  ni. i»”,-  2  2309  ni.p.h. 


5.  NUMERICAL  RELATION# 

1  lb.  «  0.453ft  kg. 
t  kg— ‘2.2040  lb. 

I  ml. ~  1,000.35  m«  8,280  ft. 
1  tn *■>:), 2808  ft. 


Best  ^a' 


\\ati\e 


Cop 


